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Abstract 

This  effort  has  focused  on  verification  and  refinement  of  the  mechanism  transitions  at 
intermediate  temperatures  in  the  disk  alloy  Rene  104,  with  the  observation  of 
microtwinning,  continuous  faulting  and  dislocation  by-pass  at  successively  higher 
temperatures  above  650°C.  Evidence  for  the  twin  initiation  process  has  also  been 
obtained  via  TEM  studies  of  specimens  interrupted  after  small  strain  levels.  A 
preliminary  model  for  the  novel  microtwinning  regime  has  been  developed  that  appears 
to  provide  reasonable  agreement  with  the  present  experimental  creep  response  for  Rene 
104  and  Rene  88.  Modeling  at  the  ab  initio ,  atomistic  and  phase  field  levels  is  providing 
important  insight  into  the  activation  parameters  associated  with  the  observed  deformation 
mechanisms,  augmented  by  insight  from  3D  atom  probe  measurements  on  y  and  y' 
composition  and  ordering.  Single  crystals  of  Rene  104  have  successfully  been  grown, 
and  a  limited  number  of  micro-tensile/compression  as  a  function  of  crystal  orientation 
whave  been  performed.  A  novel  phase  field  model  of  directional  coarsening  (rafting) 
during  high  temperature,  low  stress  creep  of  blade  alloys  has  been  developed.  This 
model  accounts  for  the  local  stress  fields  associated  with  matrix  dislocations  as  well  as 
the  lattice  misfit,  and  demonstrates  promising  qualitative  agreement  with  experiment. 

1.  Research  Objectives 

This  effort  comprises  a  coordinated  team  of  researchers  from  the  Ohio  State  University 
(under  Grant  #FA9550-05-l-0135),  University  of  Michigan  (under  Grant  #FA9550-05-l- 
0100)  and  Johns  Hopkins  University  (under  Grant  #FA9550-05-l-0102).  The  goal  of  the 
program  is  to  develop  improved  models  that  will  (a)  incorporate  more  realistic 
representation  of  the  relevant  microstructures  and  micromechanisms,  (b)  enable  modeling 
for  a  range  of  relevant  service  conditions  (c)  address  time-dependent  deformation  in  both 
disk  and  blade  alloys,  and  (d)  provide  this  information  ev  to  the  component  design 
process,  building  upon  the  paths  to  the  design  process  created  in  the  DARPA  Accelerated 
Insertion  of  Materials  (AIM)  program.  This  program  is  coordinated  with  Pro  Tresa 
Pollock  (University  of  Michigan)  is  leading  efforts  to  characterize  blade  alloy 
microstructures  and  is  providing  single  crystals  materials  for  the  rest  of  the  program  and 
Pro.  Kevin  Hemker  (Johns  Hopkins  Unviersity)  who  is  focusing  on  microtensile  tests  on 
single  crystal  variants  of  the  disk  alloy. 


2.  Accomplishments 


•  Completed  assessment  of  mechanisms  operative  during  creep  in  the  two-phase 
superalloy  R104  in  the  temperature  range  from  600-850°C.  As  a  function  of  temperature, 
there  is  a  transition  from  dislocation  activity  to  microtwinning,  extended/continuous 
stacking  faults  and  dislocation  climb  by-pass  at  the  highest  temperatures. 

•  Completed  detailed  analysis  of  microtwinning  mechanism,  which  is  the  dominant 
mechanism  at  the  upper  bound  temperature  for  service  of  disk  alloys  in  present 
applications.  Microtwining,  as  well  as  the  other  precipitate  shearing  modes,  is  controlled 
by  diffusion-mediated  reordering. 

•  Determined  the  atomic-scale  processes  and  activation  barriers  associated  with 
reordering  using  the  nudged  elastic  band  method  coupled  with  first  principles 
calculations  (VASP),  which  indicate  that  the  diffusion  coefficient  for  reordering  should 
be  similar  to  that  for  Ni-vacancy  diffusion  in  the  ordered  LI2  structure. 

•  Developed  a  phase  field  dislocation  dynamics  simulation  which  explains  the  transition 
from  dislocation  activity  to  microtwinning  as  a  function  of  precipitate  structure,  applied 
stress  and  fault  energies. 

•  Used  Z-contrast  imaging  to  detect  or  the  first  time  segregation  of  heavy  elements  to 
microtwin  and  stacking  faults  in  the  precipitates  after  creep  deformation.  The 
fundamental  role  of  these  heavy  elements  (C,  Mo  and  W)  in  modifying  these  fault 
energies  has  been  explored  using  VASP  calculations. 

3.  Program  Summary 

3.1.  Introduction 

Nickel  base  superalloys  are  an  important  class  of  materials  that  are  used  in 
demanding  high  temperature  structural  applications  such  as  the  hot  zone  of  aircraft  gas 
turbine  engines  due  to  their  highly  desirable  material  properties  such  as  high  temperature 
strength  and  their  resistance  to  creep,  thermal/mechanical  fatigue  and  oxidation.  These 
unique  properties  are  attributed  to  alloy  chemistry  and  microstructure,  which  consists 
primarily  of  ordered,  N^Al-based,  L  ^-structured,  y'  precipitates  that  are  coherently 
embedded  in  a  solid  solution  y  matrix.  Dependent  upon  how  the  material  is  processed  and 
heat-treated,  the  precipitates  size,  size  distribution  and  spacing  can  be  tailored  for  optimal 
mechanical  properties.  It  is  the  size  of  these  microstructural  features  that  will  ultimately 
control  the  microscopic  deformation  mode  and  macroscopic  deformation  response. 

High  temperature  deformation  mechanisms  of  polycrystalline  superalloys  have 
been  extensively  studied  in  the  past.  These  studies  have  uncovered  a  variety  of  different 
deformation  mechanisms  that  can  become  operative  depending  on  the  testing 
temperature,  applied  stress  and  micro  structure.  At  lower  temperatures  and  higher  stresses, 
shearing  of  the  y1  precipitates  occurs  via  the  coupled  movement  of  paired  a/2<110> 
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dislocations.  These  are  perfect  dislocations  in  the  matrix,  and  when  coupled  together  by 
an  anti-phase  boundary  (APB)  in  the  y'  phase,  they  can  also  shear  the  precipitates  without 
leaving  a  planar  fault  in  their  wake  [1-3].  A  recent  summary  of  this  mechanism  and  its 
use  in  modeling  the  yield  strength  in  polycrystalline  alloys,  in  combination  with  other 
strengthening  mechanisms  such  as  grain  boundary  and  solid  solution  strengthening,  has 
been  presented  by  Kozar,  et  al.  [4].  This  mechanism  appears  to  prevail  up  to  around 
600°C  in  polycrystalline  alloys.  At  much  higher  temperatures  (around  800°C),  climb  by¬ 
pass  of  y'  precipitates  by  individual,  unpaired  a/2<110>  dislocations  is  predominant  [5]. 
We  note  that  this  mechanism  forms  the  basis  of  several  well-established  creep  models  for 
particle  and  precipitate  strengthened  alloys  [6,  7],  With  respect  to  disk  applications; 
however,  this  climb  by-pass  regime  appears  to  be  well  beyond  the  practical  service 
temperatures. 

Several  other  mechanisms  have  been  identified  in  the  temperature  range  between 
APB  shearing  and  climb  by-pass  processes.  These  mechanisms  are  distinctive  in  that 
several  types  of  planar  defects  are  formed  in  the  y'  precipitates  and/or  the  matrix.  It  is 
emphasized  that  this  temperature  regime  appears  to  correspond  to  the  onset  of  more 
strongly  thermally-activated  deformation,  and  is  essential  to  understand  with  respect  to 
improving  the  temperature  capability  of  the  polycrystalline  superalloys.  There  are  several 
reports  of  y'  precipitates  shearing  modes  and  the  formation  of  either  superlattice  intrinsic 
stacking  faults  (SISF)  and  superlattice  extrinsic  stacking  faults  (SESF),  while  a/2<110> 
dislocations  persist  in  the  matrix  [8-11].  These  "isolated  faulting"  mechanisms  require 
dissociation  reactions  at  the  matrix  precipitate  interfaces  involving  a/2<110> 
dislocations.  Another  variant  on  the  "isolated  faulting"  mode  was  initially  documented 
by  Kear  [12]  in  which  a<  1 1 2>  perfect  matrix  dislocations  shear  the  y'  particles  by 
forming  SISFs  and  SESFs  pairs  in  the  y'  precipitates.  Extended  faults  have  also  been 
reported  either  continuously  traversing  both  matrix  and  precipitates  [8,  13,  14].  These 
extended  faults  are  typically  on  the  order  of  a  few  micrometers  in  width.  Diffraction 
contrast  analysis  indicates  that  the  faults  are  predominantly  extrinsic  in  nature  following 
tensile  creep,  although  a  significant  density  of  intrinsic  faults  are  also  observed. 

Most  remarkably,  recent  experimental  work  by  Kolbe  [15]  and  Viswanathan  et  al. 
[11,  16]  reported  on  a  rather  novel  deformation  mechanism  that  operates  during  exposure 
to  moderate  stresses  and  temperatures  in  the  range  of  650°C  to  750°C.  This  mechanism 
was  denoted  as  microtw inning  due  to  the  fact  that  it  imparts  very  thin  twins  (microtwins) 
of  the  order  on  4-50  atomic  planes  in  the  deformed  microstructure  [17].  The  microtwins 
can  be  readily  identified  with  the  help  of  conventional  and  high  resolution  TEM  imaging. 
The  microtwins  shear  both  the  y  phase  and  y'  phase  and  are  commonly  seen  to  traverse 
entire  grains.  This  process  has  also  been  observed  in  the  MC2  single  crystal  alloy  using  in 
situ  TEM  deformation  studies  at  temperature  range  800-850°C.  The  observations  indicate 
that  the  dislocations  proceed  in  viscous  glide  motion  [18].  The  microtwinning  process, 
and  its  possible  relationship  to  the  other  shearing  mechanisms,  is  the  focus  of  the  present 
paper. 


It  is  noted  here  that  the  “intermediate  temperature”  range  of  interest,  roughly 
650°C  to  850°C,  is  below  that  explored  in  many  previous  studies  of  creep  in  single 


3 


crystal  Ni-base  superalloys  for  blade  applications  [19,  20].  The  polycrystalline  Ni-base 
superalloys  considered  presently  have  been  developed  to  have  modest  lattice  misfit  and  y' 
volume  fractions  close  to  50%.  It  is  noted  that  the  mechanisms  may  differ  from  those 
observed  in  the  single  crystal  blade  alloys  where  larger  y'  volume  fractions  and  lattice 
misfit  between  the  matrix  and  precipitates  create  significant  internal  stresses  that  make 
the  interfaces  even  more  profound  barriers  to  slip.  Nevertheless,  as  discussed  below, 
important  commonality  between  mechanisms  in  disk  and  blade  alloys  may  exist  for 
similar  deformation  regimes.  Indeed  many  of  the  important  initial  studies  of  the 
“intermediate  temperature”  regime  have  been  performed  on  early  “generation”  Ni-base 
superalloys,  and  these  studies  are  discussed  in  the  context  of  more  recent  work  on 
polycrystalline  disk  alloys. 

In  section  3.2  of  this  report,  the  present  status  of  knowledge  concerning 
microtwinning  in  Ni-base  superalloys  is  reviewed  in  detail.  In  section  3.3,  new  high- 
resolution  imaging  of  the  microtwins  and  twinning  partial  dislocations  is  presented  which 
unambiguously  confirms  earlier  work  concerning  the  nature  of  the  twinning  partials. 
These  results  are  consistent  with  the  hypothesis  that  thermally-activated  microtwinning  is 
controlled  by  a  diffusion-mediated  reordering  within  the  y'  precipitates.  In  section  3.4,  a 
simple  model  of  this  process  is  presented.  The  detailed  process  by  which  reordering  may 
take  place  is  explored  using  ab  initio  calculations  in  section  3.5,  which  results  in  a 
quantitative  correlation  between  reordering  and  Ni  self-diffusion.  In  section  3.6,  evidence 
for  segregation  to  stacking  faults  and  microtwins  in  the  y'  phase  are  rpesented  and 
discussed.  Finally,  in  section  7,  the  importance  of  reordering  for  several  of  the  other 
deformation  mechanisms  observed  in  both  disk  and  blade  type  alloys  is  also  discussed. 


3.2.  Previous  Experimental  Evidence  for  Microtwinning  at  Intermediate 
Temperatures 

There  have  been  several  previous  observations  of  microtwinning  following  the 
higher  temperature  deformation  of  Ni-base  superalloys.  These  results  and  proposed 
mechanisms  are  briefly  summarized,  and  additional  evidence  provided  by  recent  high 
resolution  scanning  transmission  electron  microscopy  (HR-STEM)  observations  are 
presented. 

2.1.  Observations  of  Microtwinning  After  Creep  of  y/y'  Superalloys 

There  have  been  several  early  reports  of  the  occurrence  of  microtwinning  following  creep 
of  single  crystal  and  polycrystalline  superalloys.  Kear,  Oblak  and  Kear  [21,  22]  provided 
an  early  evidence  for  extended  stacking  faults  and  deformation  twins  in  Mar-M200  single 
crystals  oriented  near  [110]  deformation  axis.  The  formation  of  the  microtwins  was 
rationalized  in  the  context  of  y'  intrinsic  stacking  faults  (SISF)  that  are  commonly 
observed  in  the  deformed  microstructure.  It  was  postulated  that  a  build-up  of  the  faults 
on  consecutive  planes  will  lead  to  the  formation  of  twins  in  y\  As  the  built-up  becomes 
larger,  the  array  of  Shockley  partials  that  are  present  at  the  interface  will  create  a  high 
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stress  configuration.  Annihilation  of  these  partials  within  the  y  channels  will  lead  to 
formation  of  continuous  microtwins.  This  proposed  process  of  microtwin  formation  is 
illustrated  in  Fig.l. 

Similar  reports  of  microtwins  in  [1 10] -oriented  single  crystals  of  CMSX-4  after 
relatively  low  temperature  creep  (<800°C)  have  also  been  reported  by  several  other 
groups  [23,  24],  Ardakani  et  al.  [25]  provided  one  of  the  most  comprehensive  studies  of 
microtwinning  in  the  single  crystal  superalloy,  SRR-99  deformed  along  [Oil]  at  750°C 
and  850  MPa.  In  this  study,  a  range  of  twin  thicknesses  were  observed  via  TEM,  and  the 
twin  volume  fraction  was  observed  to  increase  with  increasing  strain.  Their  analysis  was 
the  first  to  show  conclusively  that  twinning  contributes  significantly  to  the  total  creep 
strain  at  smaller  strain  levels.  "Embryonic"  twins  were  proposed  to  arise  from  the 
dissociation  of  a/2[  110]  matrix  dislocations  via  the  reaction: 

a/2[l  10]  =  a/6[l  T  2]  +  a/3[121]  (1) 

This  energetically  unfavorable  dissociation  reaction  was  presumed  to  be 
stimulated  by  the  presence  of  the  y'  phase.  As  the  a/2[l  10]  matrix  dislocation  encounters 
a  y'  particle,  the  a/3[121]  dislocation  could  shear  the  particle,  forming  an  SISF  or  SESF, 
but  the  a/6[  1  12]  partial  would  be  pinned  at  the  entrance  side  of  the  particle.  This  is 

the  case  since  shearing  of  the  y'  particle  by  the  complete  a/2[  110]  displacement  would 
create  a  higher  energy  APB.  While  not  specified  in  their  paper,  the  twins  were  presumed 
to  grow  by  movement  of  successive  a/3  [121]  partials  on  adjacent  (1  1 1)  planes. 


a)  Single  fault  formation 
(planar  view) 

3B6 


b)  Overlapping  faults  (twins) 


c)  Continuous  faults  (twins) 


Fig.l.  Proposed  mechanism  by  Kear  et  al.  [21,22]  for  y'  precipitate  shear,  (a)  A 
perfect  dislocation  a/2[  112]  dissociates  into  partials  at  the  y/y'  interface;  only  the 
l/3[  1 12]  shears  the  y'  phase,  (b)  Overlapping  faults  are  formed  and  (c)  annihilation 
of  the  dislocations  at  the  y/y'  interface,  formation  of  continuous  twins. 

Microtwins  were  also  reported  in  early  work  by  Guimier  and  Strudel  [26]  on  the 
polycrystalline  Waspalloy  with  about  25  volume  percent  y'  particles.  When  this  same 
alloy  was  given  a  solutionizing  heat  treatment  and  quench  in  order  to  suppress  y' 


precipitation,  twins  were  not  observed.  This  supports  the  concept  that  twinning  is  related 
to  the  presence  of  the  y'  precipitates,  possibly  promoting  the  dissociation  reaction  (Eqn. 
1)  above. 

3.2.2.  Chen  and  Knowles  Mechanism 

Knowles  and  Chen  [27]  report  that  SESF  configurations  in  the  larger  y'  particles  coexist 
with  microtwinning  for  particular  crystal  orientations  in  CMSX-4.  The  initiation  and 
propagation  of  the  microtwins,  as  well  as  the  formation  of  isolated  SESFs  in  the  y' 
particles,  are  presumed  to  be  associated  with  the  operation  of  two  different  a/3 <11 2> 
partials  on  every  {111}  plane  (e.g.  2aB  then  2aC  as  shown  in  Fig.2a).  Passage  of  the 
first  net  2Da  will  create  an  SESF,  while  subsequent  passages  will  thicken  the  twin.  For 
an  opposite  sense  of  applied  shear  stress,  an  SISF  would  be  expected.  The  authors  make 
the  very  interesting  observation  that  microtwins  are  observed  only  in  orientations  for 
which  SESFs  are  also  observed. 


AAA 


C  C  C 
B  B  B 


B  B  B 


Fig.  2(a)  Model  of  twin  formation  by  passage  of  1/3<112>  Shockley  partials 
(2aB  and  2aC)  on  adjacent  planes.  After  Knowles  and  Chen  [9].  The  four  a/3 [112] 
partials  create  a  true  twin  in  LI 2  structure  of  y’-phase.  (b)  Model  of  twin  formation  by 
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passage  of  1/6<112>  Shockley  partials  (Da)  on  adjacent  planes  Kolbe  [15].  Thermally 
activated  reordering  is  required  to  create  a  true  twin  in  LI 2  structure  of  y’-phase. 

The  authors  propose  that  the  a/3 <11 2>  twinning  partials  arise  from  a  pole 
mechanism,  which  originates  with  the  interaction  of  two  identical  a/2<110>  matrix 
dislocations  (DC),  e.g.: 

a/2[l  10]  +  a/2[l  10]  =  a[l  10]  =  a/3  [21  l]  +  a/3[121]  (2) 

The  temperature  and  stress  dependence  of  the  process  is  associated  with  the 
critical  size  for  nucleation  of  the  twin  from  the  pole.  Since  these  are  twinning  partials  of 
the  ordered  LI2  structure,  the  subsequent  propagation  of  the  a/3<112>  partials  is 
presumably  rapid  and  not  rate  controlling  since  the  formation  of  SESF’s  and  microtwins 
by  this  mechanism  does  not  require  local  diffusion  or  atomic  rearrangement.  A  difficulty 
with  this  model  for  twinning  is  that  it  requires  successive  shear  of  {111}  planes  in 
directions  that  are  necessarily  of  small  Schmid  factor  (presuming  that  the  net  shear 
direction,  2Da,  has  large  Schmid  factor). 

In  this  and  the  earlier  models  of  twinning  involving  motion  of  a/3<l  12>  partials,  the  LI2 
lattice  of  the  precipitates  would  be  preserved  and  "true-twins"  would  be  formed  in  the 
precipitates  and  the  matrix.  In  the  next  section,  we  will  discuss  a  novel  concept  that  was 
introduced  by  Kolbe  [28]  proposing  a  quite  different  mechanism  for  the  formation  of 
deformation  twins  in  the  superalloys. 

3.2.3.  Kolbe  Mechanism  for  Microtwinning  and  the  Reordering  Hypothesis 

Kolbe  [28]  performed  in  situ  TEM  deformation  experiments  at  elevated  temperature  that 
indicated  a  distinct  transition  in  deformation  mechanism  with  temperature.  Below  about 
760°C,  dislocation  motion  was  observed  to  be  jerky,  in  which  encounters  with  y' 
precipitates  resulted  in  distinct  waiting  times.  The  shearing  dislocations  were  reported  to 
be  a/2<l  10>  type,  although  no  direct  evidence  was  provided  to  support  this  claim.  In  situ 
deformation  at  higher  temperature  gave  rise  to  a  distinctly  different  mode  of  shearing  in 
which  extended  faults  propagate  continuously  and  viscously  through  both  particles  and 
matrix.  These  extended  faults  are  associated  with  partials  that  move  in  a  correlated 
manner  as  pairs.  Kolbe  [15]  hypothesized  that  these  partials  may  be  a/6<l  12>  partials  of 
the  same  Burgers  vector  (e.g.  a  pair  of  Da  partials),  and  that  they  may  be  traveling  on 
parallel  {111}  planes,  as  illustrated  in  Fig.2b.  Without  detailed  confirmation  of  this 
hypothesis,  Kolbe  further  deduced  that  these  were  in  fact  microtwins,  and  that  the 
temperature  dependence  of  the  process  may  be  associated  with  reordering  that  would 
ensue  in  the  wake  of  the  twinning  a/6<l  12>  partials  as  they  traverse  the  y'  particles. 

The  model  of  Kolbe  is  thus  differentiated  from  other  studies  of  microtwinning  in 
the  important  aspect  of  the  nature  of  the  twinning  partial  dislocations.  Partials  of  the  type 
a/3 <11 2>  preserve  the  ordering  present  in  the  LI2  structure  of  the  precipitates,  and 
consecutive  passage  of  these  partials  on  adjacent  (111)  planes  would  create  true  twins  in 
the  y'  particles.  In  contrast,  while  the  a/6[l  1 2]  type  partials  proposed  by  Kolbe  are 
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twinning  partials  of  the  FCC  matrix,  passage  of  these  partials  on  consecutive  (111)  planes 
in  the  LI 2  structure  would  create  a  "pseudo-twin"  which  would  no  longer  have  the  LI2 
structure.  The  pseudo-twin  would  have  an  orthorhombic  structure  in  which  a  large 
number  of  unfavorable,  high-energy  nearest  neighbor  bonds  (i.e.  Al-Al  nearest  neighbors 
in  the  Ni^Al  structure)  would  be  created.  Without  presenting  concrete  supporting 
evidence,  Kolbe  proposed  that  this  high  energy  structure  would  not  be  created,  but 
instead  the  LI2  structure  of  the  twins  in  the  precipitates  would  be  restored  via  thermally- 
activated  re-ordering,  as  illustrated  in  Fig.3  in  which  four  partial  dislocations,  grouped  in 
pairs,  are  gliding  on  consecutive  (111)  planes.  There  are  Al-Al  nearest  neighbors  created 
in  the  wake  of  the  gliding  dislocations,  as  highlighted  with  the  bonds  indicated  in  the 
figure. 


In  support  of  this  hypothesis  is  the  observation  by  Kolbe  that  the  partials  tended  to 
travel  in  pairs.  A  single  a/6[l  12]  partial  would  create  a  complex  stacking  fault  (CSF)  in 
the  y  ’  structure.  It  is  not  possible  to  re-order  the  Al-Al  nearest  neighbor  bonds  at  a  CSF 
since  any  rearrangement  of  the  atoms  at  the  fault  would  create  unfavorable  bonds  in  the 
adjacent  (111)  planes.  However,  passage  of  two  Shockley  partials  would  create  a  two- 
layer  CSF  that  can  be  reordered  into  a  two-layer  superlattice  extrinsic  stacking  fault 
(SESF).  The  detailed  steps  involved  in  this  reordering  process  are  discussed  in  Section  5. 


Fig.3.  Schematic  of  the  key  aspects  of  the  microtwining  deformation  mechanism.  Shear 
of  y’-phase  is  accomplished  by  a  pair  of  1/6[1 12]  partial  dislocations.  The  highlighted  Al- 
Al  nearest  neighbors  can  be  eliminated  by  diffusion  mediated  reordering,  as  symbolized 
by  the  circular  arrows. 

3.3.  Detailed  Analysis  of  Microtwins 

The  hypothesis  of  Kolbe  required  validation  in  several  important  respects.  First  and 
foremost  has  been  to  distinguish  whether  the  twinning  partials  are  a/6[l  12]  and  a/3[l  12] 
dislocations.  Second,  is  to  establish  the  actual  structure  of  the  microtwins  in  the  y' 
particles.  These  analyses  have  recently  been  undertaken  by  Viswanathan,  et  al.  [11,  16] 
and  Sarosi,  et  al.  [29],  and  the  results  indicate  the  twinning  partials  all  have  the  same 
Burgers  vector  and  support  the  claim  that  the  twinning  partials  are  of  the  a/6[l  12]  type. 
This  work  has  been  performed  on  Rene  88DT  deformed  in  creep  at  (834  MPa)  and 
650°C,  and  very  similar  structures  have  also  been  found  in  the  alloy  Rene  104  deformed 
at  690  MPa  and  677°C.  Under  these  conditions,  microtwinning  is  the  predominant 
deformation  process  based  on  assessment  of  multiple  grains  in  the  TEM  specimens.  An 
example  of  microtwins  traversing  a  grain  in  Rene  104  is  shown  in  Fig.4(a).  The 
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corresponding  selected  area  diffraction  pattern  showing  twin  reflections  is  shown  in 

Fig.4(b). 

An  example  of  a  high-resolution  image  of  a  microtwin  stretching  across  a  y’- 
precipitate  is  shown  in  Fig.4(c).  In  this  particular  example  the  microtwin  is  14  layers 
thick.  It  is  found  that  a  majority  of  the  twins  in  the  y’  phase  have  an  even  number  of 
layers  and  that  their  thickness  can  vary,  with  most  being  in  the  range  of  4-50  layers  thick 
[17].  The  thickness  of  the  microtwin  is  defined  here  by  the  number  of  adjacent  (111) 
planes  that  have  been  sheared  by  partial  dislocations.  Thus,  for  example,  a  four-layer 
microtwin  is  formed  by  propagation  of  four  partial  dislocations  on  four  adjacent  planes. 

This  microscopy  and  diffraction  analysis  reveals  an  important  aspect  regarding 
the  internal  structure  of  the  microtwins,  that  is,  the  state  of  order.  The  two  important 
states  that  can  be  distinguished  are  “true-twin”  and  “pseudo-twin”.  In  the  pseudo-twin 
configuration,  the  twinned  part  of  the  crystal  is  structurally  a  different  phase  as  compared 
with  the  LI 2  structure.  The  pseudo-twin  has  an  orthorhombic  structure  that  contains  Al- 
A1  nearest  neighbors.  The  pseudo-twin  can  be  distinguished  due  to  the  presence  of 
additional  diffraction  intensities  in  the  electron  diffraction  patterns.  In  [110]  zone  axis 
diffraction  patterns  that  contains  the  twinning  plane,  the  set  of  (111)  reflections  that  do 
not  correspond  to  the  twinning  plane  can  be  used  to  make  this  distinction.  In  a  pseudo¬ 
twin  structure,  an  intensity  peak  has  to  be  present  at  Vi  of  the  g(lll)  for  two  of  three 
possible  pseudo-twin  variants  in  the  edge-on  configuration.  These  types  of  reflections 
have  not  been  observed  in  any  of  the  analyzed  diffraction  patterns,  which  suggests  that 
the  micro  twins  are  indeed  true- twins. 

It  should  be  noted  that  it  is  not  possible  to  confirm  from  a  single  HRTEM  image 
that  the  examined  microtwin  is  indeed  a  “true-twin”.  This  is  so  because  there  exists  one 
pseudo-twin  variant  with  apparent  mirror  plane  symmetry  in  the  [110]  examination 
direction.  It  is  therefore  important  to  have  the  twinning  partial  dislocations  included  in 
the  area  of  analysis,  which  unambiguously  reveals  the  type  of  variant  that  is  present. 
Detailed  discussion  on  this  distinction  can  be  found  in  Sarosi  et  al.  [29]. 


Fig.4.  (a)  General  view  of  deformed  micro  structure  revealing  microtwins  in  edge  on  and 
inclined  orientation,  (b)  Indexed  SADP  showing  the  presence  of  twin  reflections,  (c) 
Microtwin  in  the  y’-phase  that  has  a  thickness  of  14  layers.  The  insert  is  a  Fourier 
transform  confirming  the  presence  of  super-lattice  reflections  of  the  y\ 

The  diffraction  contrast  analyses  performed  previously  to  determine  the  nature  of 
the  twinning  partials  [11]  are  extremely  challenging  due  to  (a)  the  close  separation 
between  some  of  the  partial  dislocations,  (b)  the  complicating  imaging  effect  of  the  y' 
particles,  (c)  the  differential  thinning  of  matrix  and  particles,  and  (d)  the  obscuring  effect 
of  the  faults  that  are  connected  to  the  twinning  partials.  High-resolution  imaging  can 
provide  much  more  definitive  conclusions  regarding  the  nature  of  the  twinning  partials 
when  compared  to  diffraction  contrast  imaging.  An  example  of  a  microtwin  in  y’-phase 
that  contains  several  twinning  dislocations  is  shown  in  Fig.5.  These  are  HAADF  STEM 
images  obtained  using  an  FEI  Titan  80-300  microscope  with  cs-correction  on  the  electron 
probe  operated  at  300kV.  This  image  is  taken  close  to  a  grain  boundary  where  the 
microtwin  terminated.  The  HAADF  imaging  reveals  steps  due  to  the  twinning  partials 
along  the  length  of  the  twin.  It  starts  as  a  two-layer  fault  and  changes  to  four-layer 
microtwin  in  Fig.5(a),  and  then  from  4  to  6  layer  microtwin  in  Fig.5(b). 

A  better  representation  of  the  nature  of  the  microtwins  and  the  partial  dislocations 
can  be  obtained  from  color-coded  images  based  on  the  central/inversion  symmetry 
parameter  as  shown  in  Fig.5(c,d).  The  assignment  of  the  value  of  the  central  symmetry 
was  done  based  on  an  algorithm  implemented  in  freely  available  AtomEye  visualization 
software  [30].  The  center  of  symmetry  color-coding  (blue  is  the  symmetry  of  the  perfect 
crystal,  while  "warmer"  colors  represent  deviation  from  the  perfect  crystal  symmetry) 
enables  us  to  establish  that  the  separation  between  the  Shockley  partials  is  very  small, 
less  than  1.5  nm  in  one  case  and  virtually  “compact”  in  the  other  case.  The  high- 
resolution  images  also  unambiguously  confirm  that  the  operative  dislocations  are 
Shockley  a/6[  1 1 2]  partials.  The  drawing  of  the  Burgers  circuit,  encompassing  the 
twinned  region,  reveals  that  the  partial  dislocations  are  two  a/6  [112]  dislocations  viewed 
in  the  30°  orientation.  A  simulation  that  depicts  the  identical  displacement  and  which  can 
be  obtained  only  from  two  1/6[1 1 2]  dislocations  on  consecutive  planes  is  shown  in 
Fig.5(e,f).  Thus,  in  this  orientation,  the  a/6[l  12]  Shockley  partials  can  be  unambiguously 
distinguished  from  a/3 [112]  as  previously  discussed  by  Sarosi  et  al.[29]. 
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3.4.  Karthikeyan  Model  for  Microtwinning 

A  quantitative  model  for  microtwinning  based  on  Shockley  partial  shearing  process 
described  above,  coupled  with  the  necessity  for  reordering  of  the  twinned  regions  back  to 
the  LI2  structure,  has  been  proposed  by  Karthikeyan,  et  al.  [31].  The  model  is  based  upon 
the  assumption  that  the  cooperative  movement  of  the  Shockley  partial  pairs  can  occur 
only  when  the  energy  of  the  two-layer  fault  in  its  wake  is  reduced  via  reordering.  Since 
reordering  is  a  time  and  temperature  dependent  process,  the  twinning  partial  velocity  also 
has  these  attributes.  The  reordering  process  is  assumed  to  operate  within  the  larger 
secondary  particles,  while  the  tertiary  particles  are  assumed  to  be  small  enough  to  shear 
athermally.  These  conditions  therefore  set  upper  and  lower  bounds  on  the  applied  shear 
stress  range  for  which  the  model  is  expected  to  apply.  The  basic  geometry  of  the  model 
is  shown  in  Fig.  6. 
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Fig.5(a)  HAADF  image  depicting  presence  of  two  30°  1/6  [112]  dislocations  at  the 
microtwin  interface  and  causing  a  change  in  thickness  from  4  to  6  layers,  (b)  HAADF 
image  depicting  the  presence  of  two  1/6  [112]  dislocations  at  the  microtwin  interface, 
causing  a  change  in  thickness  from  2  to  4  layers.  (c,d)  Corresponding  RGB  colored- 
coded  images  based  on  central  symmetry  parameter  highlighting  the  interface  of  the 
microtwin.  (e,f)  Simulated  structures  of  the  twinned  region  with  two  1/6  [112] 
dislocations  whose  Burgers  is  30°  to  the  viewing  direction  (2x  30°  1/6  [1 12]). 

Consider  the  two-layer  fault  formed  in  the  wake  of  the  first  two  Shockley  partials. 
The  fault  energy  of  this  two-layered  pseudotwin  is  rpt.  Reordering  reduces  the  energy  of 
the  two-layered  fault.  If  the  Shockley  partials  were  to  shear  the  y’  slowly,  then  there  is 
sufficient  time  for  reordering  to  take  place  and  the  energy  penalty  is  reduced,  eventually 
to  that  of  a  true  twin,  T„.  Karthikeyan  [31]  assumed  that  the  energy  drop  takes  place 
exponentially  with  time  such  that: 

HO  =  (r,,  -  r„  )exP(-.o)  +  r„  (3) 


where  K  is  a  constant  determining  the  reordering  rate  and  t  is  the  time.  The  exponential 
decay  is  a  statement  of  the  assumption  that  the  pseudotwin  ->  true  twin  transformation  is 
a  first-order  reaction.  The  resulting  strain  rate  law  based  on  the  Orowan  equation  is  given 
by: 


r  =  Ptp -K -v  =  plp-blp—r 
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where  ptp  is  the  density  of  mobile  twinning  partials,  b,p  is  the  Burgers  vector  of  the 
twinning  partial,  Dord  is  the  diffusion  coefficient  for  ordering  and  x  is  the  short  range 
diffusion  length  (assumed  to  be  several  nearest  neighbor  distances,  or  ~2  b),  and ^2  is  the 
volume  fraction  of  the  secondary  y’  precipitates.  The  effective  stress,  reff,  in  the  presence 
of  tertiary  y'  precipitates,  is  given  by: 


V=T- 
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The  experimental  values  of  parameters  such  as  dislocation  density  ptp ,  volume 


fraction  of  the  secondary  y'  precipitate  that  are  critical  to  the  prediction  can  be  determined 
directly  from  TEM  observations.  The  magnitude  of  the  strain  rate  predicted  by  this 
model  is  in  good  agreement  with  the  observed  strain  rate  [31],  particularly  for  the  strain 
level  at  which  the  density  of  twinning  partials  has  been  measured.  This  agreement 
suggests  that  the  basic  kinetics  associated  with  the  twinning  partial  velocity  are 
reasonable. 


In  this  work,  Karthikeyan  assumed  that  Dord  could  be  approximated  by  the  self¬ 
diffusion  coefficient,  and  chose  the  mean  value  between  the  diffusion  coefficients  for  Ni 
and  A1  for  purpose  of  the  initial  model  evaluation.  The  validity  of  this  assumption,  the 
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exploration  of  the  actual  reordering  process,  and  its  relationship  to  self  diffusion  has  been 
evaluated  by  ab  initio  modeling  that  will  now  be  presented. 


1/6[112]  pairs 


Athermally 
sheared 
tertiary  y’ 


Fig.6.  Schematic  illustration  of  the  microtwining  mechanism  and  the  geometry  of  the 
model,  (a)  View  with  the  {111}  glide  planes  “edge-on”  and  (b)  plan  view  showing  the 
coupled  1/6<1 12>  partials  shearing  tertiary  and  secondary  y\  Other  quantities  are 
defined  in  the  text. 

3.5.  Reordering 

3.5.1.  Pathways  of  vacancy  reordering 

The  current  high-resolution  imaging  confirms  two  important  aspect  of  the  theory 
of  microtwinning  as  previously  developed  by  Kolbe  [15]  and  Viswanathan  et  al.  [11,  16]: 
(a)  the  operative  dislocations  are  Shockley  partials  1/6[1 12]  and  (b)  these  partials 
propagate  through  y’  in  closely-separated  pairs  on  a  consecutive  (111)  planes.  Images 
such  as  Fig.  5(a)  and  5(b)  enable  us  to  establish  that  the  separation  between  the  Shockley 
partials  is  very  small,  less  than  1.5  nm  in  one  case  and  virtually  “compact”  in  the  other 
case.  The  small  separation  suggests  that  the  energy  of  the  section  of  the  microtwin  in- 
between  the  two  partial  dislocations  must  be  significantly  higher  than  the  energy  of  the 
microtwin  on  either  side  of  the  partials.  The  energy  difference  of  the  different  sections 
provides  strong  energetic  counter-balance  for  the  elastic  repulsion  between  the  Shockley 
partials  and  thus  implies  that  a  diffusion-mediated  reordering  converts  the  high-energy 
pseudo-twin  on  both  sides  of  the  dislocations  into  true-twin  (or  at  least  a  lower  energy) 
configuration.  This  suggests  that  the  reordering  is  happening  in  the  immediate  vicinity  of 
the  dislocation  cores. 

Kolbe  [15]  showed  that  the  pseudo-twins  created  by  an  even  number  of  a/6[l  12] 
partials  could  reorder  into  a  true-twin  configuration  by  the  process  shown  in  Fig.7.  The 
scenario  considered  here  is  a  two-layer  microtwin,  equivalent  to  a  two-layer  CSF.  In  the 
vicinity  of  the  two  partial  dislocations  at  the  right  hand  side  of  the  figure,  the  pseudo-twin 
is  in  the  as  sheared  configuration  and  contains  the  Al-Al  nearest  neighbors  as 
schematically  depicted  with  the  grey  bonds  across  the  fault  planes.  The  area  away  from 
the  dislocation  at  the  left-hand  side  is  reordered  into  a  true-twin  configuration.  According 
to  Kolbe,  the  reordering  can  be  achieved  by  a  two-step  exchange  process  as  shown  in 
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Fig.7.  In  the  first  exchange,  Ni  and  A1  atoms  are  replaced  as  indicated  with  the  circular 
arrow  1.  The  second  exchange  involves  the  same  Ni  atom,  but  now  it  replaces  another  A1 
atom  as  shown  with  the  circular  arrow  2.  The  three  atoms  involved  occupy  the  same 
(1 10)  plane,  which  is  perpendicular  to  the  viewing  direction.  Depiction  of  the  described 
process  from  a  [111]  view  of  the  central  plane  is  shown  in  Fig.7(b,c).  For  clarity,  only 
one  of  the  atoms  above  and  below  the  central  plane  is  shown  in  this  projection  of  the 
structure.  The  atoms  involved  in  the  exchanges  are  enclosed  in  a  box.  It  is  important  to 
point  out  that  this  proposed  process,  as  well  as  an  alternative  path  described  in  [15],  do 
not  involve  an  exchange  between  nearest  neighbors.  As  a  consequence,  these  reordering 
paths  cannot  be  accomplished  directly  by  vacancy  movement,  and  are  therefore  highly 
unlikely. 

Examination  of  the  as-sheared  and  reordered  configurations  in  Fig.7  reveals  that 
there  is  an  additional,  more  direct  way  for  the  system  to  reorder  from  pseudo-  into  true- 
twin  configurations.  The  novel,  more  direct  reordering  can  be  understood  from  a 
recognition  that  the  central  plane  of  the  faulted  structure  (central  plane  in  the  wake  of  two 
a/6[  1 12]  partial  dislocations)  consists  of  a  row  of  anti-site  point  defects  (A1N|  and  NiAI)  as 
shown  in  Fig.8(a).  The  anti-site  columns  are  in  the  close-packed  [110]  direction  and  are 
perpendicular  to  the  a/6[l  1 2]  Burgers  vector.  To  reorder  such  a  system  thus  simply 
requires  an  exchange  between  A1  and  Ni  atoms  in  the  anti-site  column  as  shown  in 
Fig.8(b).  It  should  be  noted  that  changing  of  A1  and  Ni  in  any  other  closed  packed  [110] 
directions  would  not  accomplish  reordering  into  a  true-twin  configuration  as  these  are  not 
the  anti-site  columns. 
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A1  in  a  (111)  plane  below 
O  Al  in  a  (ill)  plane  above 


Fig. 7.  Two-layer  micro-twin  that  has  been  created  by  passage  of  two  1/6  [1  12]  partial 
dislocations  on  a  consecutive  planes,  as  viewed  along  the  [110].  (b)  Identical  part  of  the 
central  plane  after  being  reordered  to  SESF.  The  atom  color-coding  is  Ni-red,  Al-blue. 
(c).  The  central  plane  in  the  two-layer  microtwin.  Relative  positions  of  Al  atoms  in  the 
upper  and  lower  adjacent  plane  are  shown. 


Under  the  premise  that  reordering  of  Ni  and  Al  atoms  in  a  microtwin  is  a 
vacancy-mediated  diffusion  process,  one  can  easily  envision  several  ways  in  which  the 
vacancy  movement  can  accomplish  the  swap  between  the  Ni  and  Al  sites.  For  example, 
in  the  simplest,  most  direct  and  highly  correlated  case,  the  vacancy  would  move  in  the 
[110]  anti-site  column  as  shown  Fig.8(c).  As  the  vacancy  moves  along,  the  Al  and  Ni 
atoms  in  the  column  will  be  placed  into  the  neighboring  sites  and  thus  assume  positions 
that  are  consistent  with  a  true-twin  configuration.  There  are  many  more  viable,  but  less 
correlated,  vacancy  paths  that  will  accomplish  the  exchange  of  the  Ni  and  Al  sites.  One 
such  example  is  shown  in  Fig.8(d).  In  this  less  correlated  movement,  the  vacancy  will 
enter  the  mixed  column  via  a  NiA  anti-site,  will  accomplish  one  exchange  between  Al 
and  Ni  site  and  will  exit  from  the  mixed  column. 
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Fig.8.  (a)  Schematic  of  the  atom  arrangement  in  the  central  lane  of  a  two-layer  complex 
fault,  (b)  Direct  exchange  of  Al  and  Ni  atoms  that  allow  transformation  of  the  “as- 
sheared”  into  a  “true-twin”  configuration,  (b)  An  example  of  a  linear  movement  of 
vacancy  that  enables  an  exchange  between  Al-Ni  sites  in  an  entire  [1 10]  column,  (c)  An 
example  of  less  correlated  case  that  allow  single  exchange  between  Al  and  Ni  sites.  The 
atom  color-coding  is  Ni-red,  Al-blue.  The  grey  box  represents  a  vacancy  position  and  the 
arrows  indicate  vacancy  steps  during  reordering. 

3.5.2.  Ab  initio  calculations 

Understanding  the  diffusion-mediated  reordering  mechanism  requires  quantitative 
knowledge  of  the  vacancy  formation  AGf  and  vacancy  migration  AGm  energies  associated 
with  a  movement  along  the  viable  reordering  paths.  The  values  of  AG  and  AG  for  the 
key  vacancy  exchanges  of  nickel  and  aluminum  antisites  have  been  studied  with  ab  inito 
calculations.  For  the  purpose  of  the  ab  initio  calculations  described  here,  the  studies  were 
preformed  in  the  context  of  faulted  N^Al  super-cell,  without  considering  the  stress  field 
of  the  partial  dislocations. 

The  calculations  reported  in  this  work  were  performed  with  Vienna  Ab-initio 
Simulation  Package  (VASP)  [32-35].  The  projector  augmented  wave  (PAW)  method  in 
combination  with  Generalized  Gradient  Approximation  (GGA)  for  the  exchange 
correlation  potential  was  used  [36].  Spin-polarization  was  considered  in  the  calculations. 
For  the  calculation  of  the  migration  AGm  energies,  we  used  Climbing  Image  Nudge 
Elastic  Band  (NEB)  method  as  implemented  in  VASP  by  Henkelman  [37].  Five  images 
were  used  in  the  NEB  calculations.  The  NEB  calculations  were  performed  with  K-point 
sampling  lxlxl  and  energy  cutoff  270  eV. 

The  calculations  were  performed  on  fully  periodic  N^Al  super-cells  that  have 
axes  defined  by  crystallographic  directions  3[110]nbai  x  [112]  nbai  x  7/3[l  11]  nbai-  A 
two-layer  complex  stacking  fault  was  introduced  into  the  system  by  displacing  the  upper 
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part  of  the  crystal  on  two  consecutive  planes  by  a/6[l  12].  An  example  of  a  super-cell 
with  a  vacancy  on  a  Ni  sub-lattice  site  in  the  center  of  the  complex  fault  (two-layer 
pseudo-twin)  is  shown  in  Fig.9. 


Fig.9.  The  super-cell  used  in  the  ab  initio  calculations  in  near  [110]  orientation.  The 
central  part  of  the  super-cell  contains  a  two-layer  complex  stacking  fault  (CSF).  A 
vacancy  is  present  on  a  Ni  lattice  site  in  the  center  of  the  pseudo-twin,  as  schematically 
emphasized  with  the  black  square.  The  atom  color-coding  is  Ni-red,  Al-blue. 

3.5.3.  Energetics  of  the  reordering  steps 

T  here  are  several  vacancy  pathways  that  can  accomplish  the  exchange  of  NiA  and 
A1Ni  anitisites.  The  reordering  pathways  considered  in  this  work  assume  that  in  the  start 
and  end  point  configurations  the  vacancy  is  situated  on  a  Ni  sublattice,  unless  in  the  end 
configuration  the  vacancy  is  in  the  row  of  antisites.  In  the  context  of  a  two-layer  complex 
stacking  fault,  this  means  that  the  vacancy  is  situated  either  on  a  Ni  site  at  the  fault 
interface  or  a  Ni  site  in  the  central  plane  of  the  fault.  The  assumption  for  the  vacancy 
being  located  on  the  Ni  sublattice  (VaNl)  stems  from  previous  computational  studies  that 
predict  a  significantly  higher  concentration  of  vacancies  on  Ni  sub-lattice  for  a  range  of 
Ni3.xAli+x  stoichiometries  (x~±0.1)  [38,  39]. 

The  vacancy  mediated  exchange  of  Ni M  and  A1N  can  be  considered  as  a  minimum 
three- step  jump  process.  The  first  step  always  involves  an  exchange  of  a  vacancy  with  an 
atom  in  the  anti-site  [1 10]  column.  Depending  upon  whether  the  initial  exchange  is  with 
AlNl  or  Ni41,  the  viable  pathways  can  be  divided  into  two  main  categories. 

A  schematic  representation  of  the  first  catgeory,  where  the  jump  sequence 
initiates  with  VaNl-Ni  x!  exchange,  is  shown  in  Fig.lO(a).  Step  1  is  an  exchange  between 
VaNl-NiA1,  step  2  is  an  exchange  between  VaAl-AlNl  (places  A1  on  the  proper  sub-lattice) 
and  step  3  three  is  VaNl-Ni\  which  puts  Ni  on  the  proper  sub-lattice  and  thus 
accomplishes  reordering  of  one  Al-Ni  pair.  The  actual  representation  of  the  atomic 
arrangement  and  the  possible  sites  that  can  be  involved  in  the  exchange  are  shown  in 
Fig.lO(b).  There  are  ten  available  VaNl  sites  for  the  initial  exchange  Va  -NiAl.  Three  of 
the  sites  are  crystallographically  unique,  labeled  as  fi,  I2  and  I3.  There  are  nine  Ni  sites 
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available  for  the  final  exchange,  five  of  which  are  crystallographically  unique.  For  visual 
clarity,  only  three  of  the  five  are  marked  in  Fig.lO(b)  as  Fi,  F2,  F3. 


Fig.  10(a).  Schematic  of  the  reordering  sequence  initiated  with  an  exchange  between  Vav 
and  NiA1.  (b)  The  actual  representation  of  the  atom  configuration  and  the  possible 
reordering  pathways.  The  atom  color-coding  is  Ni-red,  Al-blue. 

An  example  of  the  energy  landscape  for  one  of  the  most  favorable  pathways  is 
shown  in  Fig.ll(a).  The  individual  steps  are  I1S1F3,  as  referenced  in  Fig.lO(b).  The  first 
step  Ij  has  a  migration  barrier  of  0.93  eV.  A  similar  value  for  migration  energy  is 
obtained  for  the  initial  exchange  I2,  while  the  migration  barrier  associated  with  exchange 
I3  is  substantially  larger  at  approximately  1.2eV.  The  subsequent  exchange  Si  is  highly 
energetically  favorable;  the  system  lowers  its  energy  by  0.297  eV.  The  energetic 
favorability  of  Si  can  be  ascribed  to  the  fact  that  one  Al-Al  nearest  bond  is  eliminated  in 
the  system.  The  migration  barrier  (AGm)  associated  with  the  jump  is  relatively  low, 
having  a  value  of  0.734eV.  The  final  step  F3  represents  a  vacancy  exchange  with  a  Ni 
atom  at  the  interface.  This  exchange  increases  the  energy  of  the  system  by  approximately 
0.2  eV.  The  migration  barrier  associated  with  the  exchange  is  0.95eV.  It  is  important  to 
point  out  that  the  exchange  F2,  which  would  allow  for  linearly-correlated  diffusion  in  the 
anti-site  columns  (as  described  in  Fig.6b),  has  a  migration  barrier  of  approximately 
1.2eV,  as  seen  in  Fig.ll(b).  Therefore,  step  F2  is  unfavorable.  Moreover,  the  overall 
energy  increase  associated  with  this  exchange  is  approximately  0.4eV  higher  than  for  the 
F3  step. 

From  consideration  of  the  possible  final  steps,  the  F2  jump  has  the  highest 
migration  barrier  and  the  highest  energy  cost,  thereby  ruling  out  the  linearly-correlated 
diffusion  process. 

The  second  possibility  for  the  reordering  jump  sequence  is  initiated  with  VaNl- 
A1Ni  exchange.  The  reordering  sequence  can  involve  three  steps,  as  shown  in  Fig.l2(a), 
four  steps  as  shown  in  Fig.l2(c)  or  any  higher  number  of  steps.  The  reordering  paths 
presented  in  Fig.l2(a,c)  correspond  to  a  ring-path  due  to  the  fact  that  step  1  creates  an  A1 
antisite  which  is  eliminated  in  the  last  step.  The  actual  representation  of  the  atomic 
arrangement  and  a  three-step  path  I1S1F1  is  shown  in  Fig.l2(b).  If  step  3  would  continue 
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in  the  antisite  row  as  F2,  the  reordering  can  be  envisioned  as  a  process  that  temporarily 
creates  A1N|  anti-sites  due  to  the  first  step. 


Migration  Barrier 


Reaction  Coordinates  [A] 

Fig.  1 1(a)  Energy  landscape  associated  with  vacancy  steps  that  accomplish  a  reordering  of 
single  Al-Ni  pair.  The  process  initiates  with  an  exchange  between  Va-Ni.  (b)  Comparison 
of  the  two  possible  energy  barriers  for  the  final  exchange  (VaNl-NiAI)  involving  either 
step  F3  (vacancy  exchange  with  an  atom  at  the  interface)  or  step  F2  (linearly  correlated 
vacancy  exchange  with  a  Ni  atom  in  the  central  plane  of  the  fault). 
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Fie.  12(a).  Three-step  reordering  sequence  that  initiates  with  an  exchange  between  Va  '■ 
A1  .  (b)  The  actual  representation  of  the  atom  configuration  and  the  possible  reordering 
pathways.  The  atom  color-coding  is  Ni-red,  Al-blue.  (c)  Four-step  reordering  sequence 
that  initiates  with  an  exchange  between  VaNl-AlN  . 


Fig.  13  displays  the  energy  landscape  for  a  jump  sequence  that  corresponds  to 
steps  I1S1F1  in  Fig.l2(b).  The  initial  exchange  Ii  is  highly  energetically  unfavorable;  as  a 
result  of  this  exchange,  the  energy  of  the  system  increases  by  approximately  0.3 5eV. 
This  energy  increase  can  be  rationalized  by  the  fact  that  new  Al  1  anti-site  is  created 
which  shares  3  Al-Al  nearest  neighbors.  However,  the  migration  barrier  (AGm)  associated 
with  the  exchange  is  approximately  0.7eV.  Similar  values  of  the  energy  balance  and  AGm 
are  obtained  for  the  other  possible  initial  jumps.  The  second  step  VaNl-NiA1  (Si)  also 
increases  the  energy  of  the  system.  The  energy  of  the  system  is  increased  by  0.362  eV. 
The  corresponding  migration  barrier  is  approximately  AGm~1.2  eV.  Only  the  last  step 
lowers  the  energy  of  the  system.  The  energy  is  lowered  by  approximately  0.65  eV,  which 
indicates  that  it  is  highly  energetically  favorable.  The  corresponding  migration  barrier  is 
only  0.18eV. 
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Migration  Barrier 


REACTION  COORDINATES  [A] 


Fig.  13.  Energy  landscape  associated  with  vacancy  steps  that  accomplish  a  reordering  of 
single  Al-Ni  pair.  The  process  initiates  with  an  exchange  between  Va-Ni. 

3.5.4.  Assessment  of  the  reordering  pathways 


The  energy  landscape  calculations  enable  us  to  predict  the  most  likely  pathway 
the  vacancy  will  take  in  reordering.  Based  on  the  presented  calculations  it  becomes 
evident  that  the  reordering  path  initiated  with  the  Va^-Ni'1  exchange  will  be  favored 
over  the  path  initiated  with  VaNl-AlNl.  This  is  so  despite  the  fact  that  the  migration  barrier 
for  the  step  1  in  the  VaNl-NiAI  exchange  is  approximately  0.25eV  higher  than  the  step  1  of 
the  VaN,-AlNl  initiated  sequence. 

The  argument  for  this  conclusion  is  based  on  the  fact  that  the  reordering  sequence 
initiated  with  VaNi-NiAl  exchange  has  a  favorable  relation  between  the  forward  and 
backward  transition  probability  of  the  involved  steps.  That  is,  each  step  of  the  sequence 
has  higher  transition  probability  than  the  previous  step  backward.  This  favorable  relation 
guarantees  that  once  the  initial  exchange  has  taken  place,  the  next  reordering  step  will  be 
favored. 

In  comparison,  this  favorable  relation  is  not  found  for  the  reordering  sequence 
initiated  with  VaN,-AlNl  exchange.  The  transition  probability  for  step  2  is  very  low  as 
compared  to  the  transition  probability  of  the  backward  step  1 .  The  actual  energy  barriers 
are  1.2eV  for  the  forward  (Si)  jump  and  0.4  eV  for  the  backward  (  Ii)  jump.  The 


difference  of  0.8  eV  allows  us  to  estimate  based  on  T  =  v  •  ext 


-A 

,  kj 


that  the  ratio 


between  transition  rates  Ii/Si  is  approximately  14,000  at  a  temperature  of  700°C,  which 
makes  the  overall  reordering  sequence  very  unlikely. 
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This  analysis  indicates  that  the  reordering  path  initiated  with  the  VaNl-NiA1 
exchange  will  be  dominant.  The  choice  of  the  last  step  in  the  sequence  has  significant 
implications  for  the  reordering  process  as  it  was  shown  that  the  exchange  labeled  as  F2 
would  provide  a  continuous,  highly  correlated  reordering.  Due  to  the  fact  that  the 
migration  barrier  associated  with  this  jump  F2  is  approximately  0.2  eV  higher  than  the 
other  possible  final  steps,  it  can  be  concluded  that  the  linear,  highly  correlated  motion  of 
vacancy  will  not  be  the  relevant  mode  of  reordering. 

The  picture  of  the  reordering  process  that  emerges  from  these  preliminary 
calculations  is  the  following  (see  Fig.14).  (a)  The  vacancy  is  at  the  twin  interface  on  Ni 
site.  The  most  likely  reordering  jump  sequence  is  such  that  the  vacancy  enters  the  mixed 
row  of  the  central  plane  though  a  Ni  site,  which  enables  the  subsequent  switch  with  Al. 
Upon  completing  this  exchange,  the  vacancy  exits  from  the  mixed  column  of  the  central 
plane,  (b)  Following  this  single  pair  reordering,  it  is  likely  than  that  the  vacancy  will  stay 
longer  at  the  interface  due  to  its  lower  energy  of  formation  in  this  plane  (as  discussed  in 
the  following  section),  (c)  The  process  is  repeated  as  the  vacancy  can  participate  in  a 
new,  analogous  reordering  process.  In  the  case  that  the  reordering  takes  place  in  the 
vicinity  of  the  previously  reordered  pair,  the  energy  of  the  reordering  steps  is  even  more 
favorable,  (d)  The  vacancy  is  available  for  further  reordering. 


Fig.14.  Sequential  picture  of  the  reordering  process  as  deduced  from  the  ab  inito  analysis, 
(a)  Vacancy  (VaNl)  reorders  one  AlNl-NiAI  pair  through  the  consecutive  jumps  in  the 
antisite  column,  (b)  Vacancy  VaNl  diffuses  at  the  interface  or  its  vicinity,  (c)  Vacancy 
(VaNl)  enables  analogous  reordering  of  Al^-Ni^  pair,  (d)  Vacancy  is  available  for  further 
reordering. 

3.5.5.  Vacancy  formation  energy 
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The  vacancy  formation  energy  AG^orm  at  the  Ni  sites  that  serve  as  the  starting 
point  for  the  reordering  was  estimated  with  ab  initio  calculations.  As  previously 
discussed,  there  exists  three  crystallographically  unique  Ni  sites  from  which  the 
reordering  sequence  can  be  started.  One  site  is  in  the  central  plane  and  two  sites  are  at  the 
microtwin  interface.  All  three  different  sites  were  previously  shown  in  Fig.  10(b)  as  being 
associated  with  exchanges  li ,  h  and  I3,  respectively. 

For  an  intermetallic  compound  such  as  M3AI,  the  quantity  that  dictates  the 
concentration  of  point  defects  is  the  effective  formation  energy  AG^orm.  The  effective 
formation  energy  on  a  Ni  site  in  stoichiometric  M3AI  can  be  expressed  according  to 
equation  6  (the  derivation  of  the  equation  can  be  found  in  the  work  of  Mishin  [40,  41]). 
The  effective  formation  energy  depends  not  only  on  the  vacancy  formation  energy  but 
also  on  the  formation  energies  of  other  point  defects  such  as  antisites  as  they  can 
transform  through  appropriate  reaction: 

^Vam  =  eVaNI  +  Eo  +  -  (^N,ai  ~  EAlNI  )  (6) 

where,  evaNi  is  a  raw  formation  energy  at  the  Ni  site,  E0  is  energy  per  atom  in  a  perfect 
NfoAl,  ENi/u  is  the  raw  formation  energy  of  Ni  on  A1  site  (Ni  antisite),  en^  is  the  raw 
formation  energy  of  A1  on  Ni  site  (A1  antisite)  and  b=l/4  [40,  41]. 

The  values  of  the  formation  energy  at  all  three  sites  and  the  formation  energy  in 
the  bulk  are  shown  in  Table  1.  The  reported  values  were  obtained  from  two  different  K- 
point  sampling  and  energy  cutoffs.  Both  calculations  are  qualitatively  consistent  and  it 
can  be  concluded  that  the  formation  energy  in  the  bulk  is  higher  than  the  vacancy 
formation  energy  at  the  other  Ni  sites.  The  calculations  performed  with  higher  K  point 
sampling  and  higher  energy  cut-off  show  the  effective  formation  energy  is  higher  in  the 
bulk  by  approximately  0.05  eV. 

It  should  be  noted  that  these  values  of  the  effective  vacancy  formation  energy  are 
lower  than  those  previously  reported  in  experimental  or  theoretical  studies,  which 
indicate  a  value  of  1.6eV  [41,  42].  It  is  important  to  point  out  that  our  calculations  were 
performed  on  large  cells  with  limited  K-point  sampling  in  order  to  make  the  calculations 
feasible.  These  calculations  are  thus  not  expected  to  provide  highly  accurate  absolute 
values.  Nevertheless,  the  relative  trends  between  the  different  configurations  are 
expected  to  be  valid. 


Table.  1  Vacancy  formation  energy  in  M3AI  two-layer  complex  stacking  fault  as 
estimated  based  on  ab-initio  calculations. 


VACANCY  POSITION 

Vacancy  formation  energy  [eV] 

K  point  [lxlxl], 

Ecut=270  eV,  Fixed 

Cell 

K  point  [1x2x1], 

Ecut=35 1  eV, Fixed 

Cell 

VaNl  :Bulk 

1.345 

1.404 

VaNI  :Intrf.(l  1 1)/Ni[l  10] 

1.294 

1.341 

VaNl  :Intrf(l  1  l)/AlNi[l  1 0] 

1.314 

1.368 

VaNI  :Cntr(l  1 1)/Ni[l  10] 

1.378 

1.346 

VaNI 

:Cntr(l  1 1  )/AlNiNi  110] 

1.409 

1.500 
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3.5.6.  Kinetics  of  the  reordering 


Finally,  the  reordering  kinetics  can  be  estimated  from  the  knowledge  of  the 
activation  energy  (AQ)  of  the  relevant  diffusion  steps.  The  activation  energy  for  diffusion 
is  determined  as  the  addition  of  migration  and  formation  energies,  A Q  -  A(/orm  +  AGm'8. 
As  previously  discussed,  the  most  critical  step  of  the  reordering  sequence  is  the  first  step. 
Thus,  the  activation  energy  of  the  entire  sequence  can  be  approximated  with  the  value  of 
the  activation  of  the  first  step. 

The  present  calculations  reveal  that  the  effective  vacancy  formation  energy  AGf<  ™ 
for  the  initial  reordering  sites  is  lower  by  approximately  0.05eV  than  the  in  the  bulk.  The 
sites  will  thus  have  increased  concentration  of  vacancies  at  the  interface,  and  will 
promote  the  reordering  kinetics. 

T  he  migration  barrier  AGm'8  for  the  first  step  of  the  most  likely  exchange  in  the 
reordering  sequence  is  ~0.93eV.  This  value  is  about  0.07eV  higher  than  the  migration 
barrier  of  0.86eV  for  the  vacancy  jump  on  the  Ni  sub-lattice  in  the  “bulk”  section  of  the 
super-cell.  The  higher  value  of  the  migration  barrier  of  the  first  jump,  as  compared  to  the 
bulk,  will  have  a  retarding  effect  on  the  reordering  kinetics. 

Note  that  the  calculated  value  of  0.86eV  for  the  migration  barrier  in  the  bulk  of 
the  super-cell  is  lower  than  the  value  of  0.97  eV  reported  in  other  ab  initio  work  in  the 
literature  [43].  As  mentioned  above,  this  difference  is  due  to  limited  K-point  sampling 
and  the  relatively  low  energy  cutoff  270  eV  presently  used.  When  the  calculation  are 
performed  on  2x2x2  N^Al  cell,  with  K  -points  sampling  of  8x8x8  and  energy  cut-off  of 
351  eV,  the  migration  barrier  of  0.97eV  is  fully  consistent  with  the  previous  ab  initio 
work  in  the  literature  [43]. 

Considering  both  contributions  from  AGform  and  AGm'8,  it  is  revealed  that  the 
differences  are  nearly  mutually  offset,  which  suggests  that  the  activation  energy  for  bulk 
and  the  first  step  in  the  reordering  sequence  are  approximately  the  same.  This  is  a  very 
important  result  since  it  indicates  that  the  reordering  kinetics  can  be  approximated  using 
Ni  self-diffusion  in  N13AI. 

This  analysis  represents  a  simplified  picture  for  the  kinetics  of  reordering.  The 
simplifications  are  most  notably  that:  (a)  This  analysis  is  based  on  the  reordering  of 
single  pair  of  AlNl  and  Ni'1  antisites  in  a  row  of  antisites.  The  single  pair  reordering 
eliminates  two  Al-Al  nearest  neighbors  across  the  fault,  but  creates  one  Al-Al  nearest 
neighbor  bond  in  the  central  row,  which  is  eliminated  only  upon  subsequent  reordering  of 
the  next  antisite  pair.  The  energies  associated  with  the  vacancy  exchanges  will  be  altered 
for  a  partially-reordered  fault,  (b)  This  calculation  accounts  only  for  the  activation  energy 
(migration  and  formation),  but  not  for  correlation  effects,  (c)  This  analysis  is  based  only 
on  vacancy  exchanges  in  a  dislocation-free  fault  region,  and  does  not  consider  the 
possible  effect  of  the  stress  and  strain  field  from  the  core  of  twinning  partials.  (d)  Actual 
commercial  alloys  contains  several  alloying  elements  that  are  “y’  formers”  [44]  that 
preferentially  locate  in  the  ordered  precipitates  and  could  modify  the  formation  and 
migration  energies  that  have  been  calculated  above  for  binary  N^Al. 

3.5.7.  Reordering  of  Other  Faulted  Structures 
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In  the  context  of  microtwining,  it  was  shown  that  the  reordering  can  proceed  in 
the  portion  of  microtwin  that  was  created  by  an  even  number  of  identical  a/6<l  12>  type 
dislocations.  It  is  now  useful  to  generalize  the  reordering  condition  in  terms  of  the  (111) 
shear  displacement  in  LI 2  crystals.  To  fully  avoid  any  Al-Al  nearest  neighbors,  it  is  found 
that  reordering  can  be  accomplished  only  on  faults  that  are  formed  by  a  net  shear 
displacement  of  n/3[l  12]  or  n[l  10],  In  addition  to  the  microtwinning  mechanism 
described  above,  there  are  therefore  several  other  dislocation  mechanisms  that  meet  the 
reordering  criterion.  Table  2  lists  several  possible  dislocation  configurations  that  create 
faults  capable  of  reordering. 


Table.2.  Example  of  sheared  configurations  capable  of  reordering 


Description 

Schematic 

Two  1/6<1 12>  on  subsequent  planes 

T 

ll/6[112]  x-4  CUT) 

Il/6[112]  W 

Two  1/6<1 12>  on  separate  planes 

ll/6[112]  ~ 

1  1/6(112]  ^ 

1/6<1 12>  and  1/2<1 10>  on  subsequent  planes 

1  1/6[1 12] 

J[  1/2[110] 

Zonal  1/6<1 12>  and  1/2<1 10>  on  subsequent  planes 

X,  1/6[121]  J[l/2[110]  ^-4 

1 1/6(211] 

Two  1/2<1 10>  on  subsequent  planes 

1  1/2[110]  ^4 

1  l/2[lT0]  W 

Configurations  that  do  not  meet  the  specified  condition  for  reordering  are  also 
worthwhile  examining.  For  example,  one  can  imagine  a  situation^  in  which  a  two-layer 
fault  is  created  by  two  different  partial  dislocations,  such  as  a/6[l  12]  and  a/6[12  1].  The 
net  displacement  from  these  two  dislocations  is  a/6[211]  and  therefore  the  system  is  not 
capable  of  full  reordering  (SESF  configuration).  Nevertheless,  this  2xCSF  configuration 
can  be  partially  reordered  by  reducing  the  number  of  Al-Al  nearest  neighbors  in  the 
central  plane  of  the  fault  from  two  to  one,  which  leads  to  a  reduction  of  the  fault  energy. 

Interestingly,  a  four-layer  configuration  created  by  a  pair  of  a/6[l  12]  and 
a/6[12  1]  can  be  fully  reordered.  However,  the  complexity  of  the  dislocation 
configuration  and  the  fact  that  the  full  reordering  can  proceed  only  after  passage  of  4 
dislocations  may  prevent  such  a  shearing  configuration  and  explain  its  absence  in 
experimental  observations. 
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3.5.8.  Implications  for  Other  Mechanisms 


The  insights  discussed  above  regarding  reordering  may  also  be  key  for 
understanding  of  a  variety  of  other  planar  deformation  mechanisms  involving  the 
shearing  of  y'  precipitates  that  have  been  reported  to  operate  at  intermediate  temperatures 
in  Ni  base  superalloys.  Like  microtwinning,  these  mechanisms  also  appear  to  progress  in 
a  viscous  manner  and  be  thermally  activated  [11,  12,  45]. 

One  such  a  deformation  mechanism  involves  the  formation  of  “isolated”  SESFs. 
In  this  mechanism,  faulting  (an  SESF)  occurs  only  in  the  y'  phase  [8,  11,  16],  The  nature 
of  the  created  faults  implies  that  the  matrix  has  been  sheared  by  a  perfect  lattice 
translation,  while  the  y'  phase  has  been  sheared  by  a  pair  of  a/6 [112]  dislocations.  The 
nature  of  y'  phase  shearing  is  identical  to  microtwinning,  indicating  that  an  identical 
diffusion-mediated  reordering  is  also  the  rate-limiting  process  in  this  case. 

The  other  most  notable  deformation  mechanisms  that  do  not  rely  on  the  operation 
of  identical  a/6[l  12]  pairs,  but  nonetheless  may  require  reordering,  are:  (a)  the  SISF 
mechanism  that  operates  by  a/3 [112]  dislocations  and  (b)  the  dislocation  ribbon 
mechanism  that  operates  by  a[l  12]  dislocations.  Only  the  dislocation  ribbon  is  further 
discussed  in  detail  since  the  SISF  mechanism  can  be  considered  as  part  of  the  dislocation 
ribbon  mechanism,  corresponding  to  the  leading  segment  of  the  dislocation  ribbon. 

The  dislocation  ribbon  mechanism  operates  under  creep  deformation  conditions  at 
approximately  750°C  and  stresses  700-750MPa  during  primary  creep  of  single  crystal 
alloys  in  near  [100]  orientation  [12,  46,  47].  Such  temperature  and  stress  conditions  are 
comparable  to  conditions  that  induce  microtwinning  in  some  microstructures  of 
polycrystalline  Rene  104  alloy.  The  identification  of  the  dislocation  ribbon  and  its 
recognition  as  the  key  deformation  mechanism  during  primary  creep  was  first  made  in 
series  of  publications  by  Kear  at  al.  [12,  46,  48].  Their  study  was  performed  on  Mar- 
M200,  which  can  be  considered  a  1st  generation  superalloy.  In  subsequent  studies  on 
newer-generation  superalloys,  the  operation  of  the  dislocation  ribbon  has  been  confirmed 
and  proposed  as  the  dominant  deformation  process  for  primary  creep  of  single  crystals 
near  [100]  orientation  [47]. 

The  Burgers  vector  of  the  dislocation  ribbon  is  a[112],  which  is  a  perfect 
translation  in  the  LI2  ordered  crystal  structure.  The  propagation  of  the  dislocation  thus 
does  not  leave  planar  defects  in  the  y/y'  structure.  However,  the  a[l  12]  dislocation  is 
widely  spread  and  presents  a  configuration  that  can  be  distinguished  by  superlattice 
intrinsic/extrinsic  stacking  fault  (SISF/SESF)  ribbons.  A  schematic  of  the  dislocation 
ribbon  is  shown  in  Fig.15.  The  individual  partial  dislocation  as  well  as  the  fault 
components  of  this  configuration  can  be  expressed  as:  a/3[  112]  &  |S7S7r|  &  a/6[  112] 
&  \APB  &  a/6[  112]  &  \SESF\  &  a/3[  U2]. 
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Fig.  1 5.  Stacking  fault  ribbon  reproduced  from  the  work  of  Rae  and  Reed  [47]. 


Although  expressed  differently,  the  previous  work  by  Kear  et  al.  [48,  49] 
recognized  that  this  deformation  mechanism  would  be  a  viscous  process  limited  by 
diffusion.  As  stated  by  Kear,  Leverant  and  Oblak  [48]  “the  motion  of  the  dislocation 
ribbon  is  most  likely  controlled  by  diffusion  since  adjustments  are  required  at  the  core  of 
the  a/3[l  12]  leading  dislocations  as  it  enters  the  y'  particles”.  In  the  following  text,  we 
show  that  these  adjustments  are  analogous  to  the  reordering  process,  a  fact  that  has  not 
been  recognized  to  date  in  the  literature. 

To  understand  the  reordering  or  “adjustment”  process,  it  is  important  to 
understand  the  dislocation  configurations  at  the  atomic  level.  In  the  early  work  of  Kear  et 
al.  [49,  50],  it  has  been  shown  that  there  exist  two  different  dislocation  configurations 
with  a[112]  Burgers  vector  that  can  produce  the  characteristic  intrinsic/extrinsic  fault. 
One  of  the  dislocation  configuration,  which  was  subsequently  found  to  be  consistent  with 
the  experimental  observations  [12],  is  shown  in  Fig.16.  The  configuration  is  complex, 
composed  of  a  series  of  partial  dislocations  gliding  on  two  adjacent  planes.  (The  leading 
dislocation  is  on  the  left-hand  side  and  the  trailing  is  on  the  right-hand  side).  The  detailed 
depiction  of  the  dislocations  is  shown  in  Fig.l6(a),  while  a  compact  notation  is  shown  in 
Fig.l6(b).  In  the  compact  notation,  the  leading  a/3[l  12]  dislocation  is  represented  as 
B  8  +  B8.  Note  that  B  8  is  meant  to  imply  that  the  dislocation  is  spread  on  two  atomic 
planes  (zonal  dislocation),  and  is  composed  of  a  8C  partial  on  one  plane  and  a  8A  on  the 
adjacent  plane.  Similarly,  the  second  a/6[l  12]  partial  is  represented  in  equivalent  way  as 
B  8. 
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Fig.  16.  The  dislocation  configuration  in  the  intrinsic/extrinsic  fault  that  are  associated 
with  the  viscous  slip,  (a)  The  detail  of  the  dislocation  configuration  (b)  simplified 
notation  showing  the  dislocation  dipole  displacement  [49]. 

The  key  components  of  the  dislocation  ribbon  are  the  square  symbols  that  indicate 
“dipole  displacements.”  These  dipole  displacements  provide  the  necessary  shift  of  atoms 
in  the  central  plane  by  a/2 [110 ]  in  order  to  obtain  the  low  energy,  Al-Al  nearest- 
neighbor-free  configuration  in  the  wake  of  the  gliding  dislocations.  There  are  three 
different  dipole  displacements  required  in  the  wake  of  the  passing  dislocations.  Kear  et 
al.  [49]  postulated  that  these  dipole  displacements  are  responsible  for  the  propagation  of 
the  dislocation  by  viscous  glide. 

The  interpretation  given  by  Kear  et  al.  [49]  is  directly  analogous  to  reordering. 
However,  their  explanation  must  be  modified  with  respect  to  the  locations  of  the 
ordering/adjustments.  This  new  interpretation  is  depicted  in  Figs.  17  and  18.  Focusing  on 
the  first  segment  of  the  ribbon,  which  is  shown  in  Fig.17,  it  is  recognized  that  reordering 
into  an  SISF  can  take  place  only  after  the  passage  of  the  entire  a/3[l  12].  In  other  words, 
reordering  can  proceed  only  on  the  two-layer  CSF(8C)/APB(BA)  segment,  but  not  on  the 
two-layer  CSF(8C)/CSF(8A)  segment.  The  two-layer  CSF(8C)/CSF(8A)  could  be 
partially  reordered,  but  there  are  two  degenerate  states,  only  one  of  which  converts  to  an 
SISF  after  the  passage  of  the  subsequent  B8.  It  is  not  currently  clear  why  there  should  be 
a  bias  towards  a  configuration  that  subsequently  transforms  to  an  SISF,  and  thus  it  is 
argued  that  the  reordering  can  take  place  only  at  the  core  of  the  B8  Cpartial,  or  after  its 
passage. 


The  reordering  process  of  two  layer  CSF(8C)/APB(BA)  fault  to  an  SISF  is 
analogous  to  that  described  above  for  microtwinning.  It  requires  a  switch  of  atoms  in  the 
row  of  anti-sites.  The  row  of  anti-sites  is  along  the  AB/BA  direction,  as  shown  in 
Fig.l7(e).  Since  the  reordering  atoms  have  a  different  bonding  environment  than 
microtwining,  it  is  likely  that  activation  barriers  will  be  somewhat  different  than  those 
calculated  in  Section  5. 
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Fig.  17.  (a)  Edge-on  schematic  view  of  the  leading  a/2  [1 12]  in  the  dislocation  ribbon.  The 
dislocation  is  split  into  a  leading  1/3[1 12]  and  a  trailing  1/6[1 12],  which  are  separated  by 
an  SISF.  (b,c,d,e)  Plan- views  of  the  atom  arrangements  for  the  three  planes  schematically 
shown  in  (a).  For  clarity,  only  A1  atoms  are  shown  for  the  upper  (U)  and  the  lower  (L) 
plane  in  order  to  indicate  the  Al-Al  nearest  neighbor  violations.  The  location  where 
reordering  must  take  place  are  indicated  by  the  arrow.  The  atom  color-coding  is  Ni-red 
“smaller”,  Al-blue  “larger”. 

A  schematic  of  the  dislocation  components  in  the  second  segment  of  the 
dislocation  ribbon,  corresponding  to  a/6[l  12],  a/3[l  12]  and  the  coupling  faults,  is  shown 
in  Fig.18.  There  are  two  separate  locations  where  the  diffusion-mediated  reordering  is 
required.  One  is  at  the  core,  or  after  the  passage,  of  the  zonal  a/6[l  12]  partial,  and 
converts  the  two-layer  CSF  into  an  SESF.  The  other  is  at  the  core,  or  after  the  passage  of 
the  a/3  [1 12]  partial,  and  converts  two-layer  APB  into  perfect  crystal.  The  reordering  of  2- 
layer  CSF  to  an  SESF  has  the  same  activation  barriers  as  previously  described  in  Section 
5.4.  However,  this  is  not  the  case  for  reordering  of  two-layer  APB  into  perfect  crystal. 
The  bonding  environment  is  different  and  thus  different  activation  energies  can  be 
expected.  The  reordering  can  take  place  only  after  the  passage  of  both  a/6[l  12]  rather 
than  after  the  first  a/6[l  12],  as  specified  by  Kear  et  al.[49]. 

The  reordering  concept  provides  a  clear  physical  picture  for  understanding  the 
dislocation  ribbon  movement.  In  this  context,  it  becomes  clear  why  all  three  dislocation 
segments  have  difficulty  in  propagating  through  the  y' -phase,  a  fact  that  is  not  well 
understood  in  the  literature  [47].  All  three  reordering  processes  proceed  in  a  different 
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bonding  environment  and  thus  it  can  be  expected  that  they  all  will  have  different  values 
for  migration  energy,  and  thus  kinetics.  It  can  be  expected  that  the  slowest  reordering  will 
dictate  the  motion  of  the  ribbon.  Further  exploration  of  the  energetics  and  kinetics  for  the 
two  additional  complex  fault  configurations  identified  for  the  stacking  fault  ribbon  (in 
addition  to  the  two-layer  CSF  explored  in  Section  5)  will  provide  valuable  insights  into 
this  shearing  process. 
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Fig. 18.  (a)  (a)  Edge-on  schematic  view  of  the  trailing  a/2  [1 12]  in  the  dislocation  ribbon. 
The  dislocation  is  split  into  the  leading  1/6[1 12]  and  trailing  1/3[112],  both  of  which 
dislocations  are  separated  by  SESF.  (b,c,d,e)  Plan-views  of  the  atom  arrangement  for  the 
three  planes  schematically  shown  in  (a).  For  clarity,  only  A1  atoms  are  shown  for  the 
upper  (U)  and  the  lower  (L)  plane  in  order  to  indicate  the  Al-Al  nearest  neighbor 
violations.  The  two  locations  where  reordering  must  take  place  are  indicated  by  the 
arrow.  The  atom  color-coding  is  Ni-red  “smaller”,  Al-blue  “larger”. 


The  operation  of  an  individual  a/3[l  12]  super-Shockley  partial  dislocation  is 
another  important  deformation  mode  at  intermediate  temperatures  [27,  45,  51].  The 
individual  a/3  [112]  dislocation  can  be  considered  as  a  subset  of  the  dislocation  ribbon, 
which  implies  that  reordering  is  also  relevant  for  the  creation  of  the  SISF.  This  fact  is  not 
recognized  in  the  literature.  It  is  conventionally  assumed  that  the  SISF  is  created  without 
the  need  for  reordering  [27],  This  would  indeed  be  the  case  if  a/3[l  12]  dislocation  moves 
in  a  compact  core  configuration  or  a  single  plane  dissociated  configuration.  However, 
there  exists  both  theoretical  and  experimental  evidence  indicating  that  a/3  [112]  partial 
undergoes  non-planar  dissociation.  Sun  et  al.  [52]  performed  high  resolution  TEM 
imaging  which  indicates  a  non-planar,  spread  dissociation.  The  non-planar  dissociation  of 
the  a/3[l  12]  partial  dislocation  has  been  also  discussed  by  Link  [51].  An  extrinsic  fault  at 
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the  head  of  the  SISF  was  interpreted  as  the  indication  of  non-planar  spreading.  The  core 
structure  of  the  a/3[l  12]  dislocation  spread  on  two  planes  bounded  by  an  SISF  is  also 
conceptually  similar  to  what  has  been  theoretically  predicted  by  Yamaguchi  et  al.  [53]. 

The  non-planar  spreading  of  the  a/3  [112]  has  been  previously  considered  as  a 
locked  configuration.  Due  to  the  original  work  by  Giamei  at  al.  [54],  it  is  customary  to 
refer  to  the  configuration  as  a  "Giamei  lock".  The  locking  indeed  takes  place  at  cryogenic 
temperature  as  discussed  by  Sun  et  al.  [52]  in  the  study  of  N^Ga  and  N^Al.  However,  at 
higher  temperature  and  low  strain  rates,  reordering  becomes  possible  and  may  enable  the 
glide  of  the  spread  configuration.  Studies  that  establish  the  detailed  aspect  of  the 
dislocation  configuration  for  intermediate  temperature  shearing  by  a/3[l  12]  partials  in  Ni 
based  superalloy  still  needs  to  be  undertaken. 

3.6.  Segregation  at  the  Microtwin  Interface 

3.6.1.  HAADF  STEM  observations 

The  HAADF  STEM  imaging  reveals  that  reordering  is  not  the  only  diffusion 
phenomenon  associated  with  the  formation  of  the  microtwins.  As  shown  in  Fig.l9(a),  the 
microtwin  interfacial  region  is  often  brighter  in  comparison  to  the  bulk  of  the 
twin/surrounding  matrix,  which  suggests  segregation  heavy  elements.  The  fact  that  the 
enhanced  intensity  is  likely  to  represent  compositional  variations  from  heavy  elements  is 
due  to  the  nature  HAADF  imaging,  where  the  intensity  can  be  taken  as  proportional  to 
the  Z'1'7'2.  The  intensity  profile  in  Fig.l9(b)  shows  that  the  intensity  enhancement  is 
detected  in  the  interfacial  plane  and  the  adjacent  (111)  planes.  The  intensity  profile  was 
obtained  as  an  average  over  2  nm  width  of  the  profile  box,  as  highlighted  in  Fig.l9(a). 


Fig.  19.  HAADF  STEM  image  depicting  a  microtwin  in  y/y'  interface.  The  image  intensity 
is  higher  at  the  microtwin  interface  in  y'  suggesting  segregation  of  heavier  elements,  (b) 
Intensity  profile  obtained  from  a  2  nm  wide  region  perpendicular  to  the  interface. 
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The  enhanced  intensity  has  been  detected  on  many  but  not  all  the  imaged 
microtwins.  This  indicates  that  the  strength  of  the  segregation  can  vary  and  is  likely 
linked  to  the  diffusion  exposure  time;  early  forming  microtwins  would  show  higher 
intensity  enhancement.  The  enhancement  of  the  interfacial  intensity  is  observed  only  in 
the  y'  phase  portion  of  the  microtwins,  not  in  the  y  matrix  phase. 

Detailed  examination  of  the  interfacial  planes  with  the  enhanced  intensity  often 
discloses  the  occurrence  of  periodic  intensity  variation  along  the  length  of  the  interface. 
Characteristic  double  period  is  detected  in  the  planes.  This  can  be  well  seen  in  the 
HAADF  image  shown  in  Fig.20(a),  which  depicts  a  presence  of  a  two-layer  fault.  The 
periodic  intensity  variation  could  be  best  confirmed  from  a  Fourier  spectrum  of  the 
intensity  profile,  as  shown  in  Fig.20(b).  An  extra  peak  at  half  reciprocal  distance  of  the 
fundamental  peak  is  indicative  of  the  intensity  variation.  The  FFT  spectrum  from  the 
central  (111)  plane  of  the  two-layer  fault  discloses  a  strong  half  peak.  This  is  however  not 
the  case  for  the  plane  at  the  fault  interface.  The  FFT  spectrum  from  the  adjacent  planes 
on  either  side  of  the  interface  also  shows  presence  of  half  period  peak.  The  strength  of 
the  peak  is  smaller  as  compared  to  the  central  plane  of  the  fault.  Examination  of  the 
planes  that  are  further  away  of  from  the  interface  does  not  reveal  the  presence  of  any 
similar  periodic  variation. 

The  intensity  variation  along  the  length  of  the  interface  suggests  that  the 
segregation  is  limited  to  one  type  of  N^Al  sub-lattice  site.  The  subtle  superlattice  contrast 
from  (100)  planes  in  the  acquired  images  can  be  used  as  reference  frame  for  the 
identification  of  the  type  of  the  sub-lattice.  However,  due  to  the  complex  composition  of 
the  y’phase  in  Rene  104,  the  scattering  strength  of  the  Ni  (100)  sub-lattice  planes  is  not 
much  stronger  than  the  mixed  Al-Ni  (100)  sub-lattice  planes,  and  therefore  the 
identification  of  the  reference  frame  is  very  difficult  in  the  as  acquired  image.  An 
example  of  a  microtwinned  region  where  the  distinction  of  the  (100)  planes  can  be 
recognized  more  clearly  is  shown  in  Fig.21.  The  twinned  region  is  the  in  the  lower  part  of 
the  image.  The  intensity  modulation  in  the  y’  matrix  in  the  upper  part  of  the  image  is  out 
of  step  with  the  intensity  variation  in  the  adjacent  plane  of  the  twin  interface.  This 
indicates  that  the  heavy  elements  occupy  the  A1  sites  at  the  plane  adjacent  to  the 
microtwin  interface.  More  detailed  analysis  of  the  analysis  will  be  presented  elsewhere. 
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Fig.20(a)  HAADF  STEM  image  of  a  two-layer  microtwin,  (b)  FFT  spectra  from  intensity 
profiles  taken  along  the  length  of  the  microtwin  interface.  Several  profiles  starting  with 
the  central  plane  and  through  the  third  plane  from  the  interface  were  tested. 
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Fig.21.  HAADF  STEM  image  of  microtwin  interface.  The  lower  portion  of  the  image  is 
fourteen  layer  microtwin,  the  upper  portion  is  surrounding  y’  phase.  The  planes  adjacent 
to  the  microtwin  interface  display  super-lattice  periodicity  which  is  out  of  phase  with  the 
periodicity  in  the  surrounding  y’,  indicating  that  heavy  elements  are  occupying  A1  sites  in 
the  plane  adjacent  to  the  microtwin  interface. 


3.6.2.  Site  preference  of  elements  in  the  y’-phase 


The  presently  available  HAADF  imaging  provides  a  qualitative  sense  for  the 
segregation  of  heavy  elements  to  A1  sub-lattice  at  the  adjacent  planes  of  the  twin 
interface,  although  the  identity  of  the  elements  has  not  yet  been  determined.  Since  the 
currently  investigated  alloy  Rene  104  is  a  multi-elemental  system,  this  could  mean 
several  different  scenarios.  Possible  segregation  scenarios  consistent  with  the  current 
observations  are  discussed  in  the  context  of  the  y'  -phase  composition  and  the  site- 
occupancy  of  the  involved  elements.  Further  experimental  work,  such  as  Atom  Probe 
Tomography,  EDX  or  EELS,  which  is  currently  not  available,  will  need  to  substantiate 
the  segregation  scenario  discussed  below. 


The  nominal  composition  of  the  alloy  and  the  composition  of  y'  and  y  are  shown 
in  Table  3  [44,  55].  For  the  purpose  of  segregation,  it  is  important  to  focus  on  the 
elemental  composition  of  y'  phase  and  the  site  preference  for  the  substitution  elements, 
including  Co,  Cr,  Ti,  Ta,  Nb,  Mo,  W,  which  are  distributed  between  the  two  available 
sub-lattice  sites  (A1  or  Ni).  The  site  occupancy  of  the  elements  in  y'  phase  has  been 
previously  studied  with  experimental  and  theoretical  approaches.  The  experimental 
measurement  based  on  atom  probe  analysis  has  shown  that  Ti,  Ta,  W,  Mo  and  Nb 
substitute  for  Al  whereas  Co  substitutes  for  Ni  [56,  57].  The  Cr  can  be  shared  between 
both  sites  in  order  to  accommodate  the  stoichiometric  balance  [56].  The  atom  probe 
results  are  fully  consistent  with  ab  initio  calculations  for  a  substitution  of  single  element 
into  a  binary  N^Al  system.  It  was  previously  shown  that  Ti  [58,  59],  Ta  [59,  60],  Mo 
[58,  61],  Nb  [58]  all  prefer  the  Al  sub-lattice  while  Co  strongly  prefers  the  Ni  site  [58, 
59].  There  are  contradictory  reports  concerning  Cr  site  preference.  The  work  by  Jiang  et 
al.  [39]  has  shown  that  Cr  prefers  the  Al  site  while  the  work  of  Suiter  [59]  has  shown  that 
Cr  has  virtually  no  site  preference. 
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Table.  3.  Nominal  composition  of  Rene  104  [54],  The  compositions  of  y  and  y'  (only 
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atoms  in  the  phase  (y  or  y'). 


The  above  analysis  suggests  that  the  HAADF  experimental  observations  can  be 
accounted  for  by  segregation  of  heavy  elements  such  as  Ta,  Mo,  Nb,  W  from  their  Al 
sites  in  the  bulk  to  the  to  Al  sites  in  the  interfacial  region  of  the  microtwins.  Nevertheless, 
such  a  scenario  needs  to  be  substantiated  by  other  experimental  measurement,  as  it  is  also 
possible  that  a  more  complex  rearrangement  is  taking  place  including  chemical  changes 
on  the  Ni  sub  lattice-site. 

We  have  performed  preliminary  ab  initio  calculations  that  show  that  the  elements 
such  as  Mo  and  Ta  have  a  significant  preference  for  the  Al  site  in  the  central  and  adjacent 
planes  of  the  two  layer  twin  (SESF).  Interestingly,  the  opposite  is  true  for  two-layer 
complex  fault.  This  suggests  that  such  a  segregation  scenario  could  take  place  only  after 
reordering.  Detailed  exploration  of  the  segregation  phenomenon  is  currently  in  progress. 


3.6.2.  Effect  of  segregation  on  dislocation  mobility 

The  mobility  of  the  twinning  partials  could  be  strongly  affected  by  the  presence  of 
the  heavy  elements  at  the  interface.  One  very  obvious  effect  is  through  the  strain  field  of 
the  heavy  atoms.  The  segregation  appears  to  be  enhanced  on  the  {111}  plane  adjacent  to 
the  twin  interface,  which  is  the  precise  atomic  plane  on  which  the  next  twinning  partial 
(of  the  subsequent  pair)  must  move.  Therefore,  the  observed  segregation  places  heavy 
elements  in  the  exact  location  to  maximize  their  benefit  for  retarding  the  motion  of  the 
partials. 


Another  effect  may  be  through  the  energy  of  the  fault  itself.  The  glide  of  two 
1/6[1 12]  partial  dislocations  on  a  pre-existing,  reordered  microtwin  with  segregation  will 
cause  the  creation  of  CSF-like  environment.  This  two-layer  CSF-like  structure  is 
energetically  less  favorable  in  the  presence  of  segregation.  Thus,  it  can  be  expected  that 
the  higher  fault  energy  can  oppose  the  movement  of  the  partials,  decreasing  the  mobility 
of  the  dislocations.  Finally,  since  the  heavy  elements  are  typically  slow  diffusers,  the 
reordering  reaction  will  possibly  be  slower.  All  of  these  potential  effects  must  be 
explored  in  more  detail  in  order  to  develop  a  more  quantitative  model  of  the  influence  of 
these  alloying  additions. 
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3.7.  Influence  of  Microstructure  and  Applied  Stress  on  Propensity  for 
Microtwinning 

The  scale  of  the  precipitate  structure  and  the  relative  volume  fraction  and  size  of 
the  precipitates  can  be  varied  via  cooling  rate  from  above  the  y'  solvus  temperature  (for 
supersolvus  heat  treated  materials),  as  well  as  by  the  subsequent  aging  treatment.  The 
terminology  adopted  in  the  literature,  the  coarser  precipitates  are  termed  "secondary" 
while  the  finest  population  are  termed  "tertiary". 


Time, hr 


Fig.  22.  Creep  data  for  Rene  88  in  two  initial  microstructural  conditions  (cooled  at  two 
different  rates  from  the  supersolvus  temperature)  showing  the  strong  effect  of  tertiary 
volume  fraction  of  creep  strength.  Coarsening  of  tertiary  population  occurs  at  these 
temperatures. 

3.7.1.  Effect  of  Microstructure  (Cooling  Rate)  on  Creep  Response 

As  exemplified  in  Fig.22,  the  scale  of  the  micro  structure  can  have  a  significant  effect  on 
the  creep  rate.  The  corresponding  energy-filtered  TEM  images  [11]  of  the  secondary  and 
tertiary  y'  are  shown  in  the  insets.  Shown  in  this  figure  are  creep  tests  performed  on  Rene 
88DT  at  838  MPa  and  650°C  for  samples  cooled  at  two  different  rates  (secant  cooling 
rates  of  205  and  24°C/s).  The  significant  effect  of  microstructure  on  the  creep  response 
is  apparent.  Clearly,  the  shape  of  the  curves,  the  overall  strain  rate  and  time  to  a  given 
strain  is  strongly  dependent  on  the  initial  micro  structure.  However,  even  more  important 
is  the  recognition  that  the  dominant  deformation  processes  are  also  distinctly  different  in 
these  specimens:  isolated  faulting  of  the  secondary  particles  for  the  coarse  micro  structure 
and  microtwinning  for  the  fine  microstructure,  as  shown  in  Fig.23.  The  former 
deformation  mode  is  characterized  by  a/2[110]  dislocation  activity  in  the  matrix,  with 
apparent  shearing  of  the  secondary  y1  precipitates  by  the  presence  of  SESFs  constrained  to 
the  particles  (Fig.23(a)).  This  deformation  mode  is  contrasted  with  that  in  the  finer,  fast- 
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cooled  microstructure  where  the  microtwinning  process  discussed  above  is  predominant. 
In  addition  to  this  microstructure  effect,  there  is  also  an  effect  of  applied  stress  on 
deformation  mechanism.  For  the  fine  microstructure,  at  higher  stress  (Fig.  23c),  a/2[l  10] 
dislocation  loops  are  observed  wrapping  many  of  the  secondary  y'  precipitates,  indicating 
that  Orowan  looping  is  prevalent.  Remarkably,  few  if  any  microtwins  are  observed,  as 
they  are  in  this  same  micro  structure  at  lower  stress  (Fig.  23b). 


Fig.23:  Different  deformation  mechanisms  in  crept  R88DT  deformed  at  650°C.  838  MPa 
(in  a  and  b)  and  976  MPa  (in  c).  (a)  In  the  coarse  microstructure  at  838  MPa,  isolated 
faulting  is  observed  in  which  1/2<[1 1 0]  dislocations  are  observed  in  the  matrix  while 
precipitates  are  sheared  leaving  SESFs.  (b)  In  the  fine  micro  structure  at  838  MPa, 
microtwinning  is  dominant  after  2%  creep  strain,  (c)  In  the  fine  microstructure  at  976 
MPa,  Orowan  looping  is  observed 

Insight  into  the  microstructural  factors  affecting  the  transition  from  matrix  dislocation 
activity  to  microtwinning  have  been  obtained  in  recent  experiments  coupled  with  phase 
field  dislocation  dynamics  modeling  [62].  These  results  are  summarized  presently.  In 
order  to  explore  the  initiation  of  the  microtwinning  process  in  more  detail,  creep  tests 
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have  been  explored  at  small  strain  level.  The  substructure  of  a  Rene  104  sample 
deformed  to  0.05%  plastic  strain  following  deformation  at  677°C  and  724MPa  is  shown 
in  Fig.24.  Under  these  deformation  and  micro  structure  conditions  (fast  cooling/fine 
microstructure),  microtwinning  is  observed  to  be  the  dominant  deformation  process  at  1 
and  2%  strain.  However,  at  0.05%  plastic  strain,  a  distinctly  different  substructure  is 
observed.  The  matrix  phase  deforms  first  at  this  early  stage  of  creep  deformation.  Fig.24 
shows  a/2<110>  dislocations  being  emitted  from  an  intergranular  M23C6  carbide. 
Detailed  dislocation  analysis  results  show  that  these  dislocations  all  have  the  same 
Burgers  vector  and  that  they  have  dissociated  into  Shockley  partial  dislocations  that  have 
a  Burgers  vector  of  the  type  a/6<112>.  The  leading  partial  dislocations  appear  to  have 
"percolated"  between  the  secondary  y'  precipitates,  leaving  them  unsheared.  While  the 
process  whereby  this  structure  subsequently  evolves  into  microtwins  is  presently  not 
clear,  the  decorrelation  of  a/2[110]  matrix  dislocations  nevertheless  appears  to  be  a 
necessary  precursor  to  microtwinning. 


Fig.  24.  Post  creep  dislocation  substructure  tested  at  677°C  and  724MPa  and  interrupted 
after  a)  0.05%  strain  showing  a/2<l  10>  dislocation  dissociation  into  Shockley  partial 
dislocations  bounding  matrix  intrinsic  stacking  faults  and  b)  0.2%  strain  showing  a  high 
density  of  matrix  intrinsic  stacking  faults  resulting  from  the  propagation  of  the  leading 
Shockley  partial  dislocation. 
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3.7.2.  Modeling  of  Dislocation/Precipitate  Interactions 


There  have  been  several  simulation  studies  of  the  interaction  of  dislocations  with 
precipitates  for  the  superalloy  structures.  The  early  work  of  Mohles  and  Nembach  [63] 
considered  the  propagation  at  strongly-  and  weakly-coupled  a/2<110>  dislocations  for 
the  case  of  relatively  low  volume  fraction  y'  microstructures.  More  recent  dislocation 
dynamics  calculations  by  Rao  et  al.  [64]  further  refined  the  modes  of  interactions  and 
strength  factors  that  can  be  derived  from  precipitate  interactions.  These  simulations 
focused  on  APB  shearing  modes. 

More  recent  simulations  by  Douin,  et  al.  [65]  provided  new  insights  into  the  role 
of  the  precipitate  structure  in  promoting  dissociation  of  matrix  dislocations.  In 
simulations  similar  to  those  of  Douin,  et  al  [65],  the  interaction  of  a/2<  1 1 0>  dislocation 
with  precipitates  and  the  various  kinetic  pathways  associated  with  their  dissociation  have 
been  explored  using  a  different  computational  approach  -  the  so-called  microscopic 
phase-field  (MPF)  model  [45.  The  MPF  enables  treatment  of  complex  ((  shearing 
processes  without  any  a  priori  assumptions  about  dislocation  geometry,  core  structure, 
and  formation  of  stacking  faults  {Shen,  2007  #42].  The  model  is  an  extension  of  the 
Peierls  model  of  dislocation  [66]  with  the  incorporation  of  the  generalized  stacking  fault 
(GSF)  energy  [67]  and  3D  anisotropic  microelasticity  theory  [68].  It  utilizes  a  description 
of  dislocations  in  terms  of  shear  type  inelastic  strain  fields,  which,  combined  with  the 
voxel-to-voxel  elastic  interaction  kernel,  provides  an  efficient  treatment  for  arbitrary 
dislocation  geometries.  The  model  stems  from  the  phase  field  dislocation  model  [60]  and 
its  modification  [69]  with  a  focus  at  the  atomic  length  scales.  With  ab  initio  GSF  energy, 
the  MPF  model  is  capable  of  treating  dislocation  core  structures  quantitatively  .  In 
addition,  various  planar  faults  are  treated  essentially  as  extended  core  structures  and  their 
formations  are  treated  in  a  variational  way.  As  applied  to  the  study  of  dislocation  activity 
in  the  y/y'  microstructure,  the  MPF  model  is  advantageous  in  handling  both  arbitrary 
configurations  of  dislocations  and  microstructure,  dislocation  curvatures,  dislocation 
dissociation  and  stacking  fault  formation,  anisotropic  elastic  interactions  among 
dislocations  and  between  dislocations  and  microstructure,  and  taking  into  account  the 
more  physical  diffuse  dislocation  cores  instead  of  singular  (Volterra  type)  cores.  As  a 
result  of  the  MPF  studies,  an  understanding  of  the  critical  microstructure  features  and 
deformation  conditions  that  favor  microtwinning  versus  the  various  stacking  faulting 
modes  will  have  important  implications  for  the  resultant  creep  strength. 

A  microscopic  phase  field  simulation  was  developed  to  study  the  experimentally 
observed  a/2<l  1 0>  dislocation  dissociation  into  Shockley  partials  since  it  is  believed  that 
this  process  is  a  precursor  to  the  microtwinning  mechanism  as  described  in  a  preceding 
section.  The  simulation  is  based  a  single  screw-oriented  dislocation  as  it  interacts  with  a 
secondary  y'  precipitate  population  with  a  size  (diameter)  of  300  nm  and  an  interparticle 
spacing  of  75  nm.  A  shear  stress  of  115MPa  acting  on  the  dislocations  in  a  direction 
parallel  to  the  burgers  vector  of  the  leading  Shockley  partial  AS  (inset  of  Fig.25)  was 
used.  The  shear  direction  maximizes  the  force  on  the  leading  partial.  Three  configurations 
were  initially  chosen  to  show  the  feasibility  of  this  model  to  capture  the  dissociation 
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process  and  the  results  are  presented  in  Fig.25(a-c)  for  the  same  total  simulation  time. 
The  first  two  cases  where  the  intrinsic  stacking  fault  energy  was  fixed  at  10mJ/m2  for 
cases  with  and  without  tertiary  y'  precipitates  are  presented  in  Fig.25(a-b). 


(c) 


Figure  25.  Phase  field  simulation  of  a/2<110>  dislocation  (AB)  interaction  with  a 
secondary  y'  precipitate  at  an  applied  shear  stress  of  115MPa,  300nm  y'  precipitate  size, 
and  75nm  y  channel  width  spacing.  The  simulation  shows  the  effect  of  tertiary  y' 
precipitates  and  stacking  fault  energy  a)  lOmJ/m  1SF  energy  without  tertiary  y'  b) 
10mJ/m  ISF  energy  w/tertiary  y'  and  c)  lOOmJ/m  ISF  energy  w/tertiary  y'.  Periodic 
boundary  condition  is  applied  in  the  vertical  direction. 


In  the  phase  field  model,  starting  with  the  case  without  tertiary  y'  precipitates, 
dissociation  into  Shockley  partial  dislocations  occurs  when  the  incoming  dislocations 
come  in  contact  with  the  narrow  y  channel.  As  the  applied  stress  is  below  the  Orowan 
stress  for  the  perfect  dislocation  (AB),  it  is  sufficient  to  drive  the  leading  partial  (A8)  into 
the  channel  because  of  the  smaller  line  tension  (which  is  approximately  proportional  to 
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b2,  where  b  is  the  Burgers  vector  of  the  partial).  Note  the  stress  resolved  on  the  trailing 
partial  is  only  half  that  on  the  leading  partial  for  the  orientation  chosen.  The  leading 
Shockley  sweeps  out  and  creates  an  intrinsic  stacking  fault  in  the  matrix  while  the  trailing 
partial  remains  at  the  entrance  to  the  y  channel.  Similar  results  are  observed  when  the 
simulation  is  run  with  a  random  distribution  of  tertiary  y'  precipitates  (Fig.  5b),  but  in  this 
case  the  presence  of  the  tertiary  y'  offers  an  additional  impediment  to  dislocation  motion 
as  APB’s  are  created  in  the  wake  of  the  passing  a/2<110>  dislocation,  and  CSF’s  are 
created  in  the  tertiary  y'  after  being  sheared  by  the  Shockley  partial  dislocations.  For  the 
same  microstructure  with  tertiary  y',  but  increasing  the  matrix  stacking  fault  energy  10 
fold  to  lOOmJ/nr  leads,  to  a  quite  different  result,  as  shown  in  Fig.25(c).  Under  these 
circumstances,  when  the  dislocation  makes  its  way  to  the  confines  of  the  y  channel,  due 
to  the  increase  in  matrix  stacking  fault  energy,  the  dislocation  does  not  dissociate  into 
Shockley  partials.  Instead,  the  dislocation  remains  intact  as  a  full  a/2<110>  and  is 
immobilized. 

These  preliminary  results  of  the  phase  field  model  demonstrate  several  important  factors 
that  trigger  the  decorrelation  of  the  a/2[l  10]  dislocations,  and  leading  ultimately  to 
deformation  via  microtwinning.  Several  conclusions  may  be  drawn  at  this  point,  and  are 
summarized  below: 

(a)  Not  surprisingly,  the  matrix  stacking  fault  energy  is  an  important  factor  determining 
how  the  secondary  y’  particles  are  by-passed,  and  dictating  the  critical  stress  for  by-pass. 
High  fault  energy  inhibits  decorrelation,  and  ultimately  microtwinning. 

(b)  At  a  stress  below  the  Orowan  bowing  stress  for  a/2[l  10]  dislocations  (appropriately 
modified  for  frictional  stresses  or  the  presence  of  tertiary  particles),  decorrelation  will  be 
favored  for  microstructures  in  which  the  spacing  between  secondary  particles  is  smaller. 
This  rationale  qualitatively  explains  the  mechanism  change  for  coarse  and  fine 
microstructures  shown  between  Figs.23(a)  and  23(b). 

(c)  If  the  stress  is  in  excess  of  the  Orowan  bowing  stress,  then  the  matrix  dislocations 
would  remain  undissociated,  regardless  of  stacking  fault  energy.  This  argument  is 
supported  by  the  mechanism  change  noted  in  Fig.23(c).  Note  also  that  the  absence  of 
twins  at  this  higher  stress  level  would  be  difficult  to  understand  if  dislocation  activity  and 
microtwinning  were  operative  in  parallel.  The  phase  field  simulations  instead  suggest 
that  if  the  microstructure  and  stress  conditions  are  favorable  to  a/2[110]  matrix  activity, 
rather  than  decorrelation,  then  microtwinning  will  not  be  operative. 

(d)  The  orientation  of  the  applied  shear  stress  relative  to  the  operative  Burgers  vector 
should  have  a  significant  role  on  decorrelation.  In  the  simulations  shown  in  Fig.25,  the 
force  acting  on  the  leading  partial  is  twice  that  acting  on  the  trailing  partial.  If  instead, 
the  applied  shear  stress  were  oriented  such  that  the  trailing  partial  were  more  highly 
stressed,  then  decorrelation  will  not  be  favored.  This  indicates  that  decorrelation  and 
microtwinning  may  only  be  favored  over  an  relatively  narrow  range  of  applied  shear 
stress  orientation. 
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3.7.3.  Subsequent  Formation  of  Twins 


The  experimental  evidence  and  modeling  just  described  indicate  that 
decorrelation  of  a/2<110>  dislocations  is  a  necessary  precursor  to  microtwinning.  A 
detailed  understanding  of  this  progression  remains  the  subject  of  investigation.  Several 
possibilities  are  now  briefly  mentioned  in  light  of  the  present  understanding.  A 
comprehensive  review  of  deformation  twinning  has  been  provided  by  Christian  and 
Mahajan  [70].  For  instance,  many  authors  have  proposed  a  pole  mechanism  to  describe 
twinning  in  fee  crystals  [71,  72].  In  these  models,  the  common  theme  for  the  formation  of 
a  twin  is  the  creation  of  a  pole  and  twinning  dislocation  from  a  perfect  a/2<110> 
dislocation.  The  pole  dislocation  essentially  acts  as  an  anchor  whereas  the  twinning 
dislocations  (a/6<112>  Shockley  partial  dislocation)  can  sweep  out  on  a  {111}  glide 
plane,  rotating  about  a  pole.  There  are  opposing  views  concerning  how  twins  form  via  a 
pole  mechanism.  This  topic  is  beyond  the  scope  of  this  paper  and  will  not  be  detailed 
here  [70]. 

In  the  model  proposed  by  Mahajan  and  Chin[73],  two  co-planar  a/2<110> 
dislocations  with  different  Burgers  vectors  are  envisioned  to  react  and  form  three 
identical  Shockley  partial  dislocations  on  adjacent  {111}  glide  planes,  which  can  then  act 
as  the  nucleus  of  a  three  layer  twin.  From  the  limited  experimental  evidence  of 
microtwinning  in  the  polycrystalline  y’  strengthened  superalloys  investigated  recently 
[11,  17]  all  the  a/2[l  10]  dislocations  that  are  observed  to  dissociate  into  Shockley  partial 
dislocations  (e.g.  see  Fig.  24)  have  the  same  Burgers  vector.  This  observation  suggests 
that  the  microtwinning  mechanism  in  these  alloys  differs  from  that  proposed  by  Mahajan 
and  Chin. 

In  the  double  cross  slip  model  originally  proposed  by  Pirouz  [74],  a  single  screw- 
orientated  a/2[  110]  dislocation  dissociates  into  two  Shockley  partial  dislocations.  Due  to 
the  larger  shear  stress  acting  on  the  leading  Shockley  partial  dislocation,  it  will  bow  out 
in  a  manner  reminiscent  of  a  Frank-Read  source  and  create  an  intrinsic  stacking  fault. 
When  the  leading  Shockley  partial  re-combines  with  the  static,  trailing  partial,  the 
a/2[l  10]  dislocation  that  is  once  again  created  can  cross-slip  (along  its  screw  orientation) 
by  a  single  {111}  plane,  then  re-dissociate  into  Shockley  partial  dislocations  in  the  same 
manner  as  the  dissociation  on  the  initial  {111}  plane.  If  this  operation  were  to  be 
replicated  on  multiple,  adjacent  {111}  planes,  then  Shockley  partial  dislocations 
possessing  the  same  Burgers  vector  would  now  be  operative  and  thus  a  microtwin  could 
potentially  develop.  Based  on  the  present,  limited  TEM  evidence,  it  is  not  possible  to 
definitively  validate  this  mechanism.  However,  it  is  noted  that  this  mechanism  is 
consistent  with  the  operation  of  single  a/2[l  10]  slip  systems  as  the  source  of  the  twins. 
Another,  closely  related  process  could  be  cross-slip  and  subsequent  dissociation  of 
a/2[l  10]  dislocations  onto  the  intrinsic  stacking  faults  in  the  matrix,  present  following 
decorrelation.  This  process  would  also  create  leading  Shockley  partial  pairs  of  the  same 
Burgers  vector  that  could  then  cooperatively  shear  the  larger  y'  precipitates  via  the 
reordering  process 
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3.8.  Conclusions 


A  summary  of  TEM  studies  of  polycrystalline  and  single  crystal  Ni-base 
superalloys  at  intermediate  temperatures  (~650-750°C)  has  been  presented.  There  is  a 
considerable  variety  of  deformation  mechanisms  that  have  been  reported  following 
deformation  in  this  regime.  It  is  becoming  increasingly  clear  that  the  mechanisms 
operative  are  dependent  on  details  of  microstructure,  as  well  as  the  stress  (or  strain  rate) 
and  temperature. 

While  APB  shearing  of  the  y'  precipitates  by  coupled  a/2<l  10>  dislocations  is  the 
predominant  process  at  low  temperatures,  in  this  intermediate  temperature  regime,  other 
shearing  modes  become  active,  including  the  formation  of  SESFs  and  microtwins.  These 
mechanisms  have  been  observed  at  relatively  low  stress  levels,  below  the  yield  strength, 
and  appear  to  depend  strongly  on  temperature.  These  are  particularly  unusual 
characteristics  for  deformation  twinning  in  general,  and  the  thermally-activated  nature  of 
this  mechanism  in  the  superalloys  is  linked  to  the  diffusion-mediated  reordering  of  the 
ordered  structure  within  the  precipitates  after  being  sheared  by  the  a/6<112>  twinning 
partials.  The  reordering  process  removes  wrong  nearest-neighbors  at  the  complex 
stacking  faults  created  by  pairs  of  twinning  partials  acting  on  adjacent  {111}  planes. 
Reordering  enables  pairs  of  a/6<112>  partials  to  shear  the  secondary  y'  particles  at 
modest  stress  levels-a  process  not  possible  in  the  absence  of  reordering.  Detailed 
analysis  of  the  atomic  steps  of  reordering  by  first  principles  calculations  coupled  with  the 
nudged  elastic  band  method  suggests  that  the  diffusion  coefficient  for  reordering  should 
be  similar  to  that  for  Ni  self-diffusion  in  the  bulk,  under  the  simplifying  assumption  of 
ordered  N^Al  being  representative  of  the  y'  precipitate  composition. 

An  important  conclusion  from  the  present  knowledge  of  the  mechanisms  that 
depend  on  reordering  is  that  they  may  initiate  at  temperatures  that  are  lower  than  those 
required  for  general  climb  of  dislocations.  Reordering  is  a  conservative  process 
involving  diffusion  distances  of  atomic  dimensions,  while  climb  by-pass  requires  long- 
range  diffusion  over  distances  comparable  to  the  precipitate  size.  Therefore,  these  novel 
shearing  processes  may  be  relevant  to  the  detailed  understanding  and  potential  future 
application  of  this  class  of  alloys  since  they  are  responsible  for  the  transition  from  the 
lower  temperature,  relatively  rate-independent  deformation,  to  the  higher  temperature, 
rate-dependent  processes. 

Finally,  the  observation  of  segregation  of  heavy  elements  at  the  fault  and  twin 
planes  within  the  y'  particles  indicates  that  this  feature  is  closely  related  to  the  reordering 
process.  The  elements  involved  in  the  segregation  have  not  yet  been  determined,  but 
their  presence  may  be  vitally  important  in  promoting  faulting  and  twinning  in  these 
systems,  as  well  as  moderating  the  rate  of  creep  through  potent  interaction  with  the 
twinning  partials. 
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